Study of the decay K + — ► ir + w in the momentum region 140 < < 199 MeV/c 
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Experiment E949 at Brookhaven National Laboratory has observed three new events consistent 
with the decay K + — » ir + vv in the pion momentum region 140 < P n < 199 MeV/c in an exposure 
of 1.71 x 10 12 stopped kaons with an estimated total background of 0.93 ± 0.17(stat.)to.24(syst.) 
events. This brings the total number of observed K + —> n + vv events to seven. Combining this 
observation with previous results, assuming the pion spectrum predicted by the standard model, 
results in a branching ratio of B(K + — » ■k + vv) = (l-73i^ 0B ) x 10 -10 . An interpretation of the 
results for alternative models of the decay K + — > ix + +nothing is also presented. 



I. INTRODUCTION 

This article is a detailed report of the final results 
from experiment E949 at Brookhaven National Labora- 
tory on the study of K + — > -k + vv in the pion momentum 
region 140 < < 199 MeV/c [1]. The observation 
of K + — > ir + +nothing, a charged kaon decay to a single 
charged pion and no other observable particles, is evalu- 
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ated within the framework of the standard model (SM) 
and in terms of alternative models. 

A. Interpretation of the decay K + — > 7T~ I ~+ nothing 

The only significant SM contribution to the exper- 
imental signature K + — > n + +nothing, where nothing 
represents experimentally unobservable particles, is 
K + — > tt + i>D where vv is v e v e , v^v^ or v r v r as discussed 
in Ref . [2] . The calculation of the branching ratio has un- 
dergone continual theoretical refinement and experimen- 
tal narrowing of the relevant input parameters since the 
first modern treatment of this process [3, 4]. A recent as- 
sessment of the prediction for the branching ratio of this 
process is B(K+ -> n+vu) = (0.85 ± 0.07) x 1CT 10 [5] 
where the quoted uncertainty is dominated by the uncer- 
tainty in the Cabibbo-Kobayashi-Maskawa quark-mixing 
matrix elements. This assessment included new small 
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corrections to the charm quark contributions to the SM 
branching ratio. 

There have been many alternatives to the SM inter- 
pretation of K + — » n + + nothing signature suggested over 
the years, including the following: 

1. New physical mechanisms contributing to 
K + — > 7r + w with the usual neutrino-antineutrino 
pairs. Many models incorporating new physics 
would result in a deviation from the SM prediction 
for K + — > n + vv. A summary of these through 
mid-2007 can be found in [6] . Since that time there 
have been new calculations of the branching ratio 
in the littlest Higgs model with T-parity [7, 8], 
the possible effects on the branching ratio of a 
heavy singlet up-quark [9] and a reassessment of 
the constraints of the Minimal Flavor Violation 
Model [10]. 

2. Cases in which the neutrino flavor is not conserved. 
There are examples stemming from extended Tech- 
nicolor [11], supersymmetry (SUSY) [12], and new 
effective four-fermion interactions involving neutri- 
nos [13]. Like most examples of lepton flavor vi- 
olation in kaon decay these tend to be small, but 
there are cases such as some types of R-violating 
SUSY [14], in which K + — ► ir + vv gives the limit- 
ing constraint on some of the couplings. 

3. Reactions in which a single unseen particle recoils 
against the tt + . These include species of axions [15], 
the familon [16], sgoldstinos [17], a gauge boson 
corresponding to a new U(l)' group [18, 19], and 
various light dark-matter candidates [20-22]. In 
general these models do not predict branching ra- 
tios; rather they use K + — > n + +nothing results to 
constrain parameters. 

4. Other exotic processes. These include the effects of 
"unparticles" , which can change the SM tt + energy 
spectrum as well as the branching ratio [23] . 



B. Previous results on K + — ► tt + vv below the K n 2 
peak 

A detailed discussion of the history of measurements 
of K + — > tt + vv was given in [2]. However most of these 
measurements were made in the kinematic region in 
which the 7T + is more energetic than the tt + from the 
background reaction K + — > tt + tt° (K^), dubbed the 
"7TfP(l)" region. By contrast fewer measurements have 
been made in the U m/D(2)" region in which the n + is 
less energetic than that from K W 2 (Figure 1). As will be 
discussed below, this region is experimentally more chal- 
lenging than the ■kvv{1) region for a stopped-kaon geom- 
etry principally because the ir + from K n2 decay can enter 
the m>D(2) region if it undergoes a nuclear interaction in 
the stopping target. 
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FIG. 1: Momentum spectra of charged particles from K + de- 
cay in the rest frame. The values in parentheses represent 
the branching ratios of the decay modes [24]. The hatched 
spectrum represents the 7r + spectrum from K + — > -k + vv de- 
cay assuming the V — A interaction. The densely hatched 
regions represent the -kvv{1) and / kvv{7) E949 signal regions. 



Among the examples of nvv(2) measurements was the 
first attempt to measure K + — > n + v9 in a heavy liquid 
bubble chamber experiment [25, 26] at the Argonnc Zero 
Gradient Synchrotron that was sensitive almost entirely 
to pion momenta below 200 MeV/c. This experiment 
achieved a 90% confidence level (CL) limit on the branch- 
ing ratio of 5.7 x 10~ 5 , assuming a pure vector spectrum 
for the 7T+. Limits of 3.1 x 10~ 5 and 2.3 x 10~ 5 were 
extracted under the assumptions of tensor and scalar in- 
teractions, respectively. 

Some features of the bubble chamber experiment are 
notable. The experiment relied on the positive ir + iden- 
tification by observation of the n — > \i — > e decay chain 
(tt + — > [i + v followed by /i + — ► e + v e T>^). Although no 
timing information was available, kinematic information 
(specifically the measured range of the tt + and the angle 
between the incoming K + and outgoing 7r+) was used 
to reject background due to K + decay- in- flight. Events 
were discarded that showed evidence of a 7r + -nucleus in- 
teraction in the form of a drastic change in ionization 
along the 7r + track or a kink in the tt + trajectory. Pho- 
ton detection with a stated inefficiency of 0.02 was used 
to veto 7T° decay products and provided additional back- 
ground suppression. 

There followed a series of scintillation counter ex- 
periments by a Chicago-Berkeley group that included 
a measurement in the range 142.7 MeV/c < P ff + < 
200.9 MeV/c [27]. This yielded a 90% CL upper limit 
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on the branching ratio of 9.4 x 10 -7 assuming a vector 
spectrum. Corresponding limits were also determined 
assuming a tensor spectrum, 7.7 x 10~ 7 , a scalar spec- 
trum, 1.1 x 10~ 6 and other possible shapes. In contrast to 
the bubble chamber experiment, the counter experiment 
made use of a delayed coincidence of 3.3 ns between the 
stopped K + and the outgoing track to suppress beam- 
related background including K + decay-in-flight. A her- 
metic 4-7T sr photon detector ~10 radiation lengths (r.l.) 
thick (4.3 r.l. along the incoming beam channel) achieved 
a measured inefficiency for tt° detection of< 2.2xlCP 5 at 
90% CL for identified K„ 2 decays [28] . As with the bub- 
ble chamber experiment, the n — > \x — > e chain was used 
for positive 7r + identification and the measured range of 
the 7r + provided the kinematic information used in the 
analysis. A subsequent experiment at KEK that probed 
the ~kvv(Y) region improved the detection and identifica- 
tion of the 7r — > /i — > e decay by using 500 MHz waveform 
digitizers [29]. 

The next attempt at a measurement in the 7rz^z/(2) 
momentum region emerged out of the first phase of 
the E787 experiment at Brookhaven National Labora- 
tory [30]. The E787 detector utilized and built upon 
concepts from the earlier experiments. This experiment 
obtained a 90% CL upper limit of 1.7 x 10~ 8 , assuming a 
V — A spectrum modified by a form factor obtained from 
K + — > ir°e + v data [31]. We henceforth refer to this form 
as the "standard model" interaction. E787 also obtained 
limits of 1.4 x 10 -8 and 2.2 x 10~ 8 , respectively, assum- 
ing pure tensor and scalar interactions using 7wP(2) data 
exclusively [32]. Adding m/&(l) data, E787 improved the 
limits to 1.0 x 10~ 8 and 1.8 x 10~ 8 , respectively [33]. 

The second generation of this experiment improved the 
SM limit in the txvv(2) region to 4.2 x 10~ 9 [34] and sub- 
sequently to 2.2 x 10~ 9 [35]. Assuming tensor and scalar 
interactions E787 ultimately obtained limits of 1.8 x 10~ 9 
and 2.7 x 10~ 9 , respectively [36]. 



II. THE E949 DETECTOR 
A. Detector description 

The E787 detector was upgraded in 1999-2000 to create 
the successor experiment E949 [37] . An extensive and de- 
tailed description of experiment E949 has been provided 
elsewhere [2]. In this Section we provide a summary de- 
scription of the detector and emphasize the features es- 
sential to the ttvv(2) region. 

E949 used an incident 710 MeV/c K + beam that was 
slowed and stopped in the scintillating fiber target as 
shown schematically in Figure 2. Observation of the de- 
cay K + — > tx^vv requires detection of the incoming K + 
and outgoing 7r+ in the absence of any other coincident 
activity. The charged pion was identified kinematically 
by kinetic energy (E^), momentum (P^) and range 
measurements and by observation of the tt — > [i — > e de- 
cay sequence. Since the K + — > tt + i>D branching ratio was 
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FIG. 2: Schematic side (a) and end (b) views of the upper half 
of the E949 detector. An incoming K + is shown traversing the 
beam instrumentation, stopping in the target and decaying to 
tv + tt°. The outgoing charged pion and one photon from the 
7r° —* 77 decay are illustrated. Elements of the detector are 
described in the text. 



expected to be at the 10~ 10 level, the detector was de- 
signed to have powerful 7r + identification to reject back- 
grounds from K+ -> fi+Vft {K^ 2 ), K + -> [i+v^ (K^) 
and K + — > [L+iPv^ (K^), photon detection coverage 
over 4-7r solid angle to reject and K + — > 7r + 7r°7 
(^27)1 and efficient identification of a single incoming 
K + to suppress beam-related background. 

The incoming charged-particle beam, containing ap- 
proximately three K + for every 7r + , traversed a Cercnkov 
counter, two stations of beam wire proportional chambers 
(BWPCs), a passive BeO degrader, an active degrader 
(AD) and a beam hodoscope (B4) as shown in Figure 2. 
The BWPCs, located between the UPV and BeO, are 
not explicitly shown in the figure. Typically, 1.6 x 10 6 
kaons per second entered the target during a 2.2 s spill. 
Cercnkov photons emitted by an incoming K + {tt + ) pass- 
ing through a lucite radiator were transmitted (internally 
reflected) into 14 "kaon" ("pion") photomultiplier tubes 
to form Ck (Ct) coincidences. The photomultiplier tube 
(PMT) signals were split and fed to a discriminator and 
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a x 10 amplifier. The discriminator output was used as 
input to the time-to-digital converters (TDCs) and to the 
trigger (Section II B). The amplifier outputs were fed to 
500 MHz charge-coupled devices (CCDs) [38]. The first 
(second) BWPC station was located downstream of the 
Cerenkov counter at 168.5 (68.5) cm from the target en- 
trance and each contained three planes with sense wires 
in the vertical and ±45° (±60°) to the vertical direction. 
The wire spacing in the first (second) station was 1.27 
(0.80) mm. The BWPCs enabled detection of multiple 
beam particles. The degraders were designed such that 
incident kaons stopped within the fiducial volume of the 
scintillating fiber target. The AD consisted of 40 layers of 
2-mm thick plastic scintillator (13.9 cm diameter) inter- 
leaved with 39 disks of 2. 2-mm thick copper (13.6 cm di- 
ameter). The scintillator was azimuthally divided into 12 
sectors that were coupled by wavelength-shifting (WLS) 
fibers to PMTs that were read out by analog-to-digital 
converters (ADCs), TDCs and CCDs. These devices en- 
abled measurement of activity in the AD coincident with 
the incoming beam and outgoing products of K + decays. 
The B4 hodoscope downstream of the AD had two planes 
of 16 segmented plastic scintillator counters with 7. 2-mm 
pitch oriented at ±33.5° with respect to the horizontal 
direction. The cross-section of each counter was in a "Z 
shape" to minimize inactive area traversed by the beam 
and to improve the spatial resolution [39]. Each counter 
was connected to a PMT by three WLS fibers and each 
PMT was read out by an ADC, a TDC and a CCD. The 
B4 enabled a measurement of the target entry position of 
the beam particle as well as identification of the incident 
particle by energy loss. 

The target was composed of 413 scintillating fibers 3.1- 
m long with a 5-mm square cross-section packed to form 
a 12-cm-diameter cylinder. A number of smaller (1-, 2- 
and 3. 5-mm square) "edge" fibers filled the gaps at the 
outer edge of the target. Each 5-mm fiber was connected 
to a PMT and the output PMT signal was split and in- 
put into an ADC, a TDC, and low-gain(xl) and high- 
gain(x3) CCDs. The target fiber multiplicity and energy 
sum were also generated for triggering purposes. Multi- 
ple edge fibers were ganged onto 16 PMTs with similar 
readout. Analysis of the 500 MHz sampling information 
provided by the target CCDs was essential for isolating 
and suppressing backgrounds in the 7iw(2) region. Two 
cylindrical layers of six plastic-scintillation counters sur- 
rounding the target defined the fiducial volume. The 
inner layer of counters (dubbed "I counters" or "ICs") 
were 6. 4- mm thick with an inner radius of 6.0 cm and ex- 
tended 24 cm from the upstream end of the target. The 
5-mm thick outer scintillation counters (VC) overlapped 
the downstream end of the ICs by 6 mm and extended 
196 cm further downstream. The VC served to veto par- 
ticles that exited the target downstream of the IC. Each 
IC and VC element was instrumented with a PMT and 
read out by an ADC, a TDC and a 500 MHz transient 
digitizer(TD) [40]. 

The origin of the E949 coordinate system was the cen- 



ter of the cylindrical volume defined by the ICs. This 
point also coincided with the center of the drift chamber. 
E949 employed a right-handed Cartesian coordinate sys- 
tem with +z in the incident beam direction, +y vertically 
upward and the polar angle 9 defined with respect to the 
+z axis. The entire spectrometer was surrounded by a 1 
T solenoidal magnetic field in the +z direction. 

The drift chamber, also called the "ultra thin chamber" 
(UTC) [41], was located just outside the IC, extended ra- 
dially from 7.85 cm to 43.31 cm and served to measure 
the trajectory and momentum of the charged track from 
the target to the range stack as shown in Figure 2. Each 
of the three superlayers of the UTC contained four layers 
of axial anode wires that provided xy position informa- 
tion and two cathode foil strips that provided z position 
information. Beginning at an inner radius of 45 cm, the 
range stack consisted of 19 layers of plastic scintillator 
counters and double-layer straw chambers (RSSC) [42] 
embedded after the 10 th and 14 th layers of scintillator. 
The range stack enabled the measurement of the range 
and energy of the charged particle, the observation of 
the 7r — > [i — > e decay sequence and the measurement 
of photon activity. The 19 layers of plastic scintillator 
counters were azimuthally segmented into 24 sectors as 
shown in Figure 2. Layers 2-18 (19) were 1.9 (l.O)-cm 
thick and 182 cm long and were coupled on both ends to 
PMTs through lucite light guides. The trigger counters 
(T counters) in the innermost layer served to define the 
fiducial volume for K + decay products and were 6.4-mm 
thick and 52-cm long counters coupled to PMTs on both 
ends by WLS fibers. The T counters were thinner than 
layers 2-19 to suppress rate due to photon conversions. 
Signals from each range stack PMT were passively split 
1:2:2 for ADCs, discriminators and fan-in modules. The 
discriminator outputs were sent to TDCs and used in 
the trigger. The fanned-in analog sum of four adjacent 
sectors (dubbed a range stack "hextant") was fed into 
a single TD and was provided to the trigger. The TDs 
digitized the charge in 2 ns intervals with an 8-bit resolu- 
tion. The 500 MHz sampling permitted the observation 
of a 7r + — > fi + decay with a 5-ns separation between the 
stopped pion and the emitted muon. 

Identification of K + — > tt + vv decays required detec- 
tion of any activity coincident with the charged track. 
Photons from K n2 and radiative kaon decays were de- 
tected in a hermetic photon veto system with 4-7r sr 
solid angle coverage as shown in Figure 2. Except for the 
end caps, all photon veto detectors were lead-scintillator 
sandwich-style electromagnetic calorimeters. Other de- 
tector elements, such as the range stack, target and AD, 
also served as photon veto detectors. The barrel veto 
(BV) and barrel veto liner (BVL) covered 2/3 of 4-7r sr 
outside the range stack with a radial thickness 14.3 and 
2.29 r.l., respectively. The downstream and upstream 
end caps (ECs) consisted of 13.5 r.l. thick undoped Csl 
crystals and covered approximately 1/3 of 4-n sr [43, 44]. 
The 3.1 r.l. thick upstream photon veto was mounted 
just downstream of the Cerenkov counter. The upstream 
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FIG. 3: (Upper) Contribution of each photon veto detector in 
radiation lengths as a function of cos 6 for a photon emitted 
from the origin of the detector coordinate system. (Lower) 
Contribution of each photon veto element in the upstream 
and downstream regions within 18° of the beam. 



and downstream collar (CO) counters shown in Figure 2 
provided approximately 4.5 and 9 r.l. at normal inci- 
dence, respectively. An additional collar counter (/xCO) 
was installed downstream of the downstream CO between 
the inner face of the magnet end plate and the target [2] . 
The downstream photon veto (DPV) provided 7.3 r.l. of 
coverage downstream of the target, EC and collar. The 
AD was 6.1 r.l. thick and contributed important photon 
veto coverage in the poorly instrumented region near the 
beam axis. The thickness in radiation lengths of the pho- 
ton veto system as a function of the cosine of the polar 
angle is shown in Figure 3. 



B. Trigger 

The trigger system for E949 was designed to select 
K + — > tt + vv events from the large number of K + decays 
and scattered beam particles by requirements on the 7r + 
range, evidence of a 7r + — > decay in the range stack, 

absence of other detector activity at the time of the tt + 
and the presence of a preceding K + . The elements and 
architecture of the two-stage trigger system have been 



described previously [2]; here we describe the features 
essential for the -kvv{2) region. 

The logical OR of the following two signal triggers was 
used for the irvv(2) analysis, 



TRIGtt^(I) = KB ■ (T • 2 • IC) ■ DC • (6 ct + 7 ct ) ■ 197 t 
■zjrf ■ LOrrl ■ HEX 
{BV + BVL + EC) ■ Ll.l ■ LI. 2 (1) 



TRIGtt^(2) = KB ■ (T • 2 • IC) ■ DC • 3 ct • 4 ct • 5 ct • 6 



cl 



•(13 ct + • • • 18 ct ) • 19 ct • L0rr2 ■ HEX 
{BV + BVL + EC) ■ Ll.l ■ L1.2 . (2) 

We collectively refer to the OR of the TRIG7wz?(l) and 
TRIG7r^(2) triggers as -kvv{1+2). 

The K + beam condition KB required a coincidence 
of at least five Ck PMTs, the B4 hodoscope and the 
target with at least 20 MeV of deposited energy. The 
KB signal served as the beam strobe for the trigger. 
T ■ 2 ■ IC required a coincidence of the first two range 
stack layers in the same sector with at least one IC to 
ensure that a charged track exited the target and en- 
tered the range stack. The delayed coincidence {DC) 
required the IC time to be at least 1.5 ns later than 
the Ck coincidence to select kaon decays at rest. The 
"cf designation refers to the range stack T ■ 2 sector and 
the next two adjacent sectors that would be traversed by 
a positively charged particle in the magnetic field. For 
the TY(IGitvv{2) trigger, the charged track requirements 

3 c t • 4 ct • 5 ct • 6 ct • (13 ct H 18 ct ) • 19 ct ensured hits in 

range stack layers T through 6 to suppress contributions 
from 3-body K + decays and vetoed on hits in the outer 
layers to suppress long-range charged tracks beyond the 
'kvv{2) kinematic region. The TRIGTrz/i^l) trigger con- 
dition zfrf required the z position of the charged track 
to be within the fiducial region of all traversed range 
stack layers. The LOrrl and L0rr2 were refined require- 
ments of the charged track range taking into account the 
number of target fiber hits and the track's z position in 
range stack layers 3, 11, 12, 13 as well as the deepest 
layer of penetration in order to reject long range tracks 
such as the fi + from K + — > ii + v^ decay. The BV, BVL, 
EC and HEX requirements vetoed events with photons 
in the BV, BVL, EC and range stack, respectively. The 
Ll.l used the ratio of the height and area of the pulse(s) 
recorded by the TD to select the two-pulse signature of 
the 7t + — > /i + decay in the range stack counter in which 
the charged track was determined to have stopped. The 
LI. 2 used data digitized by the range stack ADCs to re- 
ject events with hits near the stopping counter that could 
falsely satisfy the Ll.l and to reject events with hits in 
both of the two adjacent hextants when the T ■ 2 and 
stopping counter were in the same sector. For the final 
60.6% of the data taking, an online pion Cerenkov veto 
was included in the ttvv{2) trigger to mitigate the effect 
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of an increased pion flux caused by reduced electrostatic 
separator voltage [2]. 

A subset of the data selected by the 7n/P(l+2) trigger 
is shown in Figure 4. 
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FIG. 4: Range in plastic scintillator vs. momentum for 
charged particles accepted by the nvi>(l+2) trigger. The con- 
centrations of events due to the two-body decays are labeled 
A'^-2-peak and Jf M 2-peak. The decays K + — > 7r°/i + ^ M and 
K + — > /i + ^ M 7 contributed to the muon band. The pion band 
resulted from K + — > 7r + 7r°7 decays, K V 2 decays in which the 
7r + scattered in the target or range stack and beam 7r + that 
scatter in the target. The boxes at low and high momentum 
represent the signal regions for this analysis and the previous 
■kvv(X) analysis [2], respectively. This distribution represents 
0.13% of the total kaon exposure. 



In addition to TRIGttw(2) and TRIGttw(I), "moni- 
tor" triggers were formed for calibration, monitoring, and 
acceptance and background measurements [2] . The mon- 
itor triggers selected events due to K^ and K n 2 decays 
as well as scattered beam pions (7r scat ). An additional 
"CEX" monitor trigger requiring two T ■ 2 hits was used 
to collect events resulting from the charge-exchange pro- 
cess K + n — > pKg followed by Kg — > 7r + 7r~ . Information 
derived from this CEX monitor data was used as input 
to simulation to determine the background from kaon 
charge-exchange reactions as described in Section III C 4. 
In order to measure the efficiency of the T ■ 2 • IC condi- 
tion (Section III D), we also defined a KB monitor trigger 
that required the KB condition described previously. All 
monitor triggers were prescaled to reduce their contribu- 
tion to the deadtime. 



III. DATA ANALYSIS 

The total exposure for this analysis was 1.71 x 10 12 
stopped kaons corresponding to 1.43 x 10 s ttvv{1+2) 
triggers. The total exposure was slightly less than the 
1.77 x 10 12 stopped kaons used for the nvv(l) analy- 
sis [2] because some data were discarded due to more 
stringent requirements on the reliability of the BWPCs, 
the Cerenkov counter and the target CCDs. 



A. Overview 

Identification of the K + — > n + vv decay involved pos- 
itive observation of the K + and daughter tt + in the 
absence of coincident detector activity. The -KVviV) 
and nvv(2) regions in E949 extended from 211 to 229 
MeV/c [2] and 140 to 199 MeV/c in vr+ momentum be- 
low the K n 2 peak, respectively (Figure 4). 

The 7iw(2) region potentially has a larger acceptance 
than itvv{l) because the phase space is larger and the loss 
of 7r + due to nuclear interactions in the detector is smaller 
at lower pion energies. These factors partially mitigated 
the loss of acceptance due to additional requirements 
needed to suppress background in the 7r^P(2) region. 
Compared to the previous irvv(2) analyses [35, 45], the 
acceptance was increased by enlarging the size of the sig- 
nal region. 

In a further enhancement to the previous ixvv(2) analy- 
ses [35, 45], the signal region was sub-divided into regions 
with differing signal-to-background ratios. The signal-to- 
background of each region was taken into account in the 
evaluation of B{K + —> n + vi>) using a likelihood method 
(Section HIE). 



1. Kaon-decay background 

In the 7rz/£/(l) region, the background was dominated 
by Kn2-, -K^i2, K^-y and K^ decays and was sufficiently 
suppressed by positive identification of the n + based on 
kinematic properties, observation of the 7r — » /i — > e se- 
quence and by the hermetic photon veto capability [2]. 
Previous studies [35, 45] in the nup(2) region identified 
the main background as due to K^2 decays in which the 
charged pion scattered in the target, lost energy and fell 
into the signal region. The scatter reduced the direc- 
tional correlation between the charged and neutral pions. 
Thus the photons from ir° decay were directed away from 
the high efficiency barrel region of the photon veto. This 
background was suppressed, in part, by recognition of the 
scattering process in the target. A background contribu- 
tion due to scattering of the charged pion in the range 
stack was suppressed by the track pattern and energy de- 
posit in the range stack. The photon veto served to sup- 
press these "X 7r 2-scatter" backgrounds as well as back- 
ground due to the radiative decay K^-y Background 
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due to K + — > TT + Ti~e + v e (K e ^) was suppressed by iden- 
tification of additional particles in the target. Kaon de- 
cays with a muon in the final state (K^ 2 , K^-y and K^) 
were suppressed by kinematics and the recognition of the 
7r — » /i — » e signature as well as the photon veto for the 
latter two decays. 



sample that was used to determine the selection criteria. 
We then obtained unbiased background estimates by ap- 
plying the finalized selection criteria to the remaining 
"2/3" sample of tti/D(1+2) triggers. 

The preferred method of background estimation em- 
ployed the bifurcation method illustrated in Figure 5. 
The parameter space of two sets of uncorrelated cuts 



2. Beam-related background 

The beam-related backgrounds were categorized as be- 
ing due to CEX, one beam particle (single-beam back- 
ground), and two beam particles (double-beam back- 
ground). The CEX background occurred due the produc- 
tion of a K° in the target from the charge-exchange pro- 
cess K + n —=> pK°. If the K° turned into a K\ that sub- 
sequently underwent semileptonic decay (K^ — » -K + t~v 
with £~ = e~ or /i~), the ir + could fall in the m>i>(2) 
kinematic region. CEX background was rejected by ob- 
serving the gap between the kaon and pion fibers due 
to propagation of the non-ionizing K^, by the inconsis- 
teny between the energy deposited by the K + and the 
reconstructed z of the outgoing pion and by identifica- 
tion of the accompanying negative lepton. In addition, 
requirements on the delayed coincidence between the K + 
and 7r + suppressed CEX background due to the short K\ 
flight time. 

Single-beam background was due to a K + entering 
the target and decaying in flight to produce a 7r + in the 
ixvv(2) region. Incoming beam 7r + misidentified as K + 
and scattering in the target also contributed to the single- 
beam background. Positive identification of the incoming 
particle as a kaon as well as requirements on the delayed 
coincidence between the incoming and outgoing tracks 
suppressed the single-beam background. 

The two processes (kaon decay-in-flight and pion 
scattering) that contributed to single-beam background 
formed the double-beam background when preceded by 
an additional incoming kaon whose decay products were 
undetected. Double-beam background was suppressed by 
requiring an absence of activity in the beam detectors in 
coincidence with the n + detected in the range stack. 



3. Analysis method and strategy 

We used analysis procedures and strategies similar to 
that of the E949 analysis of the ttvv(1) region [2] with 
some modification that took into account the difficulty 
of isolating some background samples in the data in the 
ixvv(2) region. As with the previous analysis, we adopted 
a "blind" analysis method in that we did not examine the 
pre-defined signal region until all background and accep- 
tance analysis was completed. Since we also attempted to 
obtain all background estimates directly from the data, 
we inverted at least one selection criteria ( "cut" ) when 
we used the ttvv(1+2) data to avoid examining the signal 
region. Every third 7iw(l+2) trigger formed the "1/3" 
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FIG. 5: Schematic of the bifurcation method. The back- 
ground level in region A can be estimated from the number of 
events observed in regions B, C and D assuming CUT1 and 
CUT2 are uncorrelated. See text for details. 

"CUT1" and "CUT2" can be divided into the four re- 
gions shown in the figure by the application of each cut 
or the inverted cut. The number of events in the sig- 
nal region "A" can be determined by application of both 
CUT1 and CUT2. If the background rejection of CUT1 
was independent of CUT2, then the number of events in 
A can be estimated as the number of events in region 
B times the ratio of the number of events in regions C 
and D or A = BC/D. In practice, we employed two 
branches for the bifurcation analysis. The "normaliza- 
tion branch" analysis was performed to obtain the num- 
ber of events, N noTm , in region B. A "rejection branch" 
analysis was used to obtain D/C. We defined the rejec- 
tion as R = (C + D)/C and obtained the background 
estimate as 



b = f x N norm /(R - 1) 



(3) 



where / = 3(3/2) for the 1/3(2/3) sample. For all 
background estimates in this analysis, the normalization 
branch was taken from the irvv(l+2) data. We used the 
irvD(l+2) data to obtain the rejection branch for all back- 
grounds except for the CEX, K e 4 and -K^2 7 backgrounds 
that could not be cleanly isolated in data. For these back- 
grounds, EGS4-based simulations [46] were employed. 
When no events (N norm = 0) were available in the nor- 
malization branch, we assigned N noTm = 1. 

We checked the validity of the background estimates 
by loosening cuts and comparing the predicted number 
of events just outside the signal region with observations 
(Section IIIC9). In addition we examined events failing 
only a single major selection criteria to search for un- 
foreseen background sources and coding mistakes (Sec- 
tion IIIB8). 
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B. Data selection 

1, Event reconstruction 

Event reconstruction was performed in a number of 
steps consisting of track-fitting in various detector sys- 
tems such as the beam-line detectors, the range stack, 
the UTC and the target. Multiple iterations of the track- 
fitting were performed in many of the detector systems 
using progressively better information from track-fitting 
from other detector systems as constraints. Events were 
reconstructed as described in [2] except as noted below. 

The following discussion focuses on the target track- 
fitting to clearly define the target-fiber classification 
scheme for use in the description of the target CCD fitter 
and the cuts that used target fiber information. In con- 
trast to the analysis of the 7ri/j/(l) region [2], the fit to the 
UTC track did not include information from the target 
fibers. Performing the target fit separately improved the 
ability to detect a pion scatter in or near the target. 

After the range stack and UTC track fitting were per- 
formed, target fibers were clustered into K + and tt + 
paths based on geometry, energy and timing informa- 
tion as shown in Figure 6. The pion fibers had to lie 
along a strip (typically 1 cm in width) along the UTC 
track extrapolated into the target, have an energy be- 
tween 0.1 and 10.0 MeV and be in coincidence with the 
reconstructed time of the tt + in the range stack (t rs ). For 
the first iteration, the kaon fibers had to have greater 
than 4 MeV of energy and be coincident with the beam 
strobe time. In subsequent iterations, fibers of lower en- 
ergy which were contiguous with the putative kaon track 
could be classified as kaon fibers. Any fiber that did not 
fall into the kaon or pion fiber categories was classified 
as a photon fiber if it had more than 0.1 MeV of energy. 
The K + decay vertex fiber was identified as the kaon fiber 
closest to the extrapolated UTC track and farthest from 
the xy position of the B4 hit. Hit fibers that were located 
on the opposite side of the decay vertex with respect to 
the outgoing track were classified as "opposite-side pion" 
fibers and were possibly due to K + decays with multiple 
charged particles or a photon conversion. The energy- 
weighted average times of the K + and 7r + hits were tx 
and t„, respectively. The sums of the K + and n + ener- 
gies were Ek and respectively. The pion and kaon 
fiber energy distributions are shown in Figure 7. 

Identified pion fibers were subjected to a least-squares 
fit to the hypothesis of a positively charged pion 
track [36]. The "target-track fitter" tested the consis- 
tency of the energy-loss-corrected UTC track with the 
energy in the fibers or the distance to the track if the 
fiber was not on the track. 

To aid in the identification of pion scattering in the tar- 
get, the activity in each of the target fibers as recorded 
by the low-gain and high-gain CCDs (Section II A) was 
fitted using a single-pulse and a double-pulse hypothesis. 
For each CCD on each fiber, a standardized K + pulse 
used for the target CCD fitter was created using kaon 
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FIG. 6: Example target xy display showing the assignment of 
pion, kaon, photon and opposite-side pion fibers as described 
in the text. The arc that terminates in the reconstructed kaon 
vertex fiber represents the extrapolated UTC track. The cir- 
cles on the track represent the UTC hits with the radius giving 
the drift distance. Only the innermost superlayer of the UTC 
is shown. The IC and VC elements with hits are indicated. 
The position of the incoming K + as reconstructed by the B4 
hodoscope information is also indicated in the figure as are the 
kaon fiber "tips" described in Section IIIB2. This event was 
selected in the K e i normalization branch (Section IIIC2). 
The measured time, energy and apparent curvature of the 
contiguous fibers classified as "photon" fibers are consistent 
with a positron and the time and energy of the opposite-side 
pion fibers are consistent with a negative pion. 



fiber data from monitor trigger data. For each fiber 
having an energy greater than a fiber-dependent thresh- 
old, typically 2 (0.5) MeV for low (high) gain, the fitting 
procedure was performed on the low-gain and the high- 
gain CCD information independently. The first step of 
the procedure was a least-squares fit to a single-pulse hy- 
pothesis for each fiber channel passing the above criteria. 
The single-pulse fit used two parameters, the pulse am- 
plitude and the time. If the probability of \ 2 CP(x 2 )) °f 
the single-pulse fit was less than 25%, a double-pulse fit 
was performed. The double-pulse fit used four parame- 
ters, the amplitudes and times for the first and second 
pulses. 



2. Requirements on it in target 

Numerous requirements were placed on the activity in 
the target to suppress background and ensure reliable 
determination of the kinematic properties of the charged 
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FIG. 7: The energy per pion fiber or kaon fiber in events in 
the if 7 r2-peak region, defined in Section IIIB4, in irvi>(l+2) 
triggers. The height of the pion energy distribution was nor- 
malized to that of the kaon energy distribution for the pur- 
pose of display. The average number of pion (kaon) fibers per 
event in the selected events is 10.2 (5.3). Note the logarithmic 
abscissa. 



pion. These requirements were based on the results of 
the target CCD fitter, the reconstructed energy and time 
of the pion and kaon fibers, the pattern of kaon and pion 
fibers relative to information from the rest of the detector 
and the results of the target-track fitter. 

Target-pulse data analysis 

Detection of pion scattering in the target in the iden- 
tified kaon fibers required reliable results from the target 
CCD fitter. An algorithm, based on the energy in the 
kaon fiber as measured by the ADC and by the time dif- 
ference t n — txi determined if the information from the 
high-gain CCD, the low-gain CCD or a combination of 
the two should be used for each fiber with an acceptable 
double-pulse fit. If V{x 2 )s of the fits for both the single- 
and double-pulse hypotheses were less than 5 x 10 -5 in 
any of the kaon fibers, then the event was rejected. In ad- 
dition, the fitted time of the first pulse {t\) was required 
to be consistent with the average time of the kaon fibers. 
If any kaon fibers failed the requirement — 6 < t\ — tx < 7 
ns, then the event was rejected. This requirement was 
made on the fitted time t\ of the single-pulse hypotheses 
if the probability of \ 2 was greater than 25% and on the 
fitted time t\ from the double-pulse hypothesis otherwise. 
The requirement on t\ rejected events in which the CCD 
fitter attempted to fit a fluctuation in the tail of the data 
pulse or when there was a large second pulse and the fit- 
ter mistakenly identified it as the first pulse. For events 
passing these criteria, the second-pulse activity in a kaon 
fiber as found by the target CCD fitting was required to 
be below 1.25 MeV when the fitted second-pulse time ti 
satisfied the coincidence condition —7.5 < t% — t n < 10 



ns. An example of the fit for a high gain CCD target 
element is shown in Figure 8. In the following we refer 
to these requirements on the CCD pulse fitting as the 
"CCDPUL" cut. 
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FIG. 8: CCD pulse fit example. The histogram represents the 
pulse height distribution for the high gain CCD data. The his- 
togram was terminated at 42 ns due to a software cutoff. The 
dashed (solid) line represents the fitted total pulse shape for 
the single-(double-)pulse hypothesis. The filled area repre- 
sents the fitted second pulse for the double-pulse hypothesis. 
The arrow indicates t ra , the expected time of the second pulse 
based on the reconstructed tt + in the range stack. This event 
was rejected because the 2.5 MeV of the fitted second pulse 
was coincident with t rs . 



Kaon fiber timing 

The target kaon fiber hits were required to be con- 
sistent with a kaon approaching the K + decay vertex. 
This consistency was enforced by requiring P(x 2 ) to be 
greater than 5% for fits to the kaon fiber hit times vs. xy 
distance to the decay vertex and vs. range. This require- 
ment removed events in which the kaon decay vertex was 
incorrectly assigned. 
Pion fiber energy 

Pion fibers were required to have energies less than 
3.0 MeV. This suppressed ir + target-scatters since the 
expected mean energy deposited in a pion fiber was ap- 
proximately 1.2 MeV. This cut had an acceptance factor 
of 89.6% (Section III D 2) due to the Landau distribution 
that describes the ionization energy deposit. 

The measured range and energy of the pion and the 
pion momentum were required to be consistent with that 
expected for a tt + using a cut on a likelihood function. 
The likelihood function was calibrated using 7r sca t mon- 
itor trigger events. In addition, the total energy of the 
pion target fibers was required to be in the range of 1 to 
28 MeV and the total energy within ±4.0 ns of t rs in the 
target edge fibers was required to be less than 5.0 MeV. 
Pattern of kaon and pion fibers 

Events with a minimum distance between the centers 
of the closest pair of kaon and pion fibers greater than 0.6 
cm, more than one fiber width, were rejected. This cut 
suppressed the CEX background. A more stringent ver- 
sion of this cut that also required that no photon fibers 
filled the gap between the kaon and pion fibers was devel- 
oped to define the normalization branches for the CEX 
(Section III C 4) and double-beam (Section III C 6) back- 
ground measurements. 



10 



Two conditions were used to enforce consistency 
among the positions of the kaon decay vertex, the kaon 
and pion clusters, and the beam particle in the B4 ho- 
doscope. The first condition required that the distance 
in the xy-plane between the hit position in the B4 ho- 
doscope and the nearest tip of the kaon fiber cluster be 
less than 1.8 cm. The kaon cluster tips were defined to 
be the two kaon fibers farthest apart from each other 
(Figure 6). The second condition required that the dis- 
tance in the xy-plane between the kaon decay vertex and 
the nearest kaon cluster tip was less than 0.7 cm. This 
requirement suppressed K^2 target-scatter background 
when the scattered n + did not emerge from the fiber 
containing the K + decay. 

The total energy of opposite-side pion fibers within 
±4.0 ns of tn was required to be less than 1.0 MeV to sup- 
press background due to K e 4 decays as well as K n 2 scat- 
ters. Hereafter, this cut is referred to as "OPSVETO" . 
Target-track fitter 

The track determined by the target-track fitter was 
required to be consistent with the information in the tar- 
get fibers and the fitted UTC track in order to suppress 
backgrounds due to pion scattering, CEX, K e 4 or a sec- 
ond beam particle in the target. For three contributions 
xl , x\ an d X? to the x 2 f° r the target-track least-squares 
fit, V{x\ + X% + Xt) was required to be greater than 1%. 
These were defined as follows: 

X% was assigned a contribution for each pion fiber tra- 
versed by the track based on the comparison of the 
observed energy with the expected energy from the 
calculated range of the track and the track momen- 
tum. 

X% was assigned a contribution based on the minimum 
distance between the track and the nearest point 
of each fiber that was traversed by the track, but 
had no observed energy. This assignment acted to 
force the fitted track to go between fibers and thus 
provided precise position information on the track. 

Xt was assigned a contribution for pion fibers that 
were not traversed by the fitted track based on the 
distance to the nearest corner. 

Events were rejected if any single pion fiber contributed 
more than 35 units to x 2 which might indicate a pion 
scatter in that fiber. The fitted target track was also 
required to intersect the kaon vertex fiber. The angle be- 
tween the reconstructed target track and the UTC track 
was required to be less than 0.01 radian at the radius 
of the IC when the range of the tt + in the target was 
less than 2.0 cm. In addition, the position of the re- 
constructed 7r + trajectories from the target and UTC 
fits were required to be well-matched at the target edge. 
Events with a kink in the target n + track were suppressed 
by requiring that the difference in the distances in the xy- 
plane of the farthest and nearest pion fibers to the center 
of the fitted helix of the UTC track was less than 0.35 
times the pion range in the target. 



3. Pion track requirements 

Good pion track reconstruction was required based on 
the x 2 °f the UTC track fit. The cut on the x 2 was depen- 
dent on the number of anode and cathode hits assigned 
to the fitted track as well as on the number of unused 
anode and cathode hits. The criteria were determined 
using both K„2 and 7T sca t monitor trigger data such that 
the 7r + momentum resolution of 2.3 McV/c for the K V 2 
peak [2] was maintained while retaining high efficiency. 

Range-stack quality cuts were placed on the probabil- 
ity of x 2 °f the range-stack track fit and the agreement 
of the z position of the extrapolated UTC track with 
the range-stack timing information and, when applica- 
ble, the RSSC information. The RSSC was not available 
for charged particles that stopped in range stack layers 6 
through 10. 




Extrapolated z at target (cm) 

FIG. 9: Extrapolated target z distribution of the charged 
track. The "7r-scatter" and "signal-like" events are taken from 
7r scat and isT M 2 monitor trigger data, respectively. The required 
minimum on the extrapolated z position and the upstream 
end of the target are indicated in the figure. 

Kaons that came to rest in the target were required 
to have energy loss in the B4 hodoscope and the tar- 
get consistent with the measured K + stopping position 
based on the z position of the extrapolated UTC track. A 
likelihood function based on these three quantities, was 
calibrated using monitor trigger data. The require- 
ment on this "Beam Likelihood" function suppressed K V 2 
target-scatter and CEX backgrounds as well as back- 
ground due to an incoming beam pion that scattered in 
the target. In addition, we required the z of the extrap- 
olated UTC track to be greater than —5.0 cm as shown 
in Figure 9. Since other cuts more effectively suppressed 
background from the downstream portion of the target 
(10 < z < 20 cm), no explicit additional cut on z was 
necessary. 
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4- Decay pion kinematic requirements 

The total range (kinetic energy) of the ir + track was 
calculated as the sum of the measured range (energy) in 
the target, IC and range stack. The total momentum 
was obtained from the curvature of the fitted track in 
the UTC corrected for energy loss in the target and IC. 
Tiny corrections were applied to R v , E n and P-^ to take 
into account the inactive material in the UTC [2]. 

The upper limit of the signal region in range, energy 
and momentum was increased with respect to the pre- 
vious wis 9 (2) analyses [35, 45] to be approximately 2.5 
standard deviations from the K^i peak similar to the 
approach used for the E949 iwv(l) analysis [2]. The 
lower limits were not changed with respect to the pre- 
vious ttvv(2) analyses. The standard signal region was 
140 < P n < 199 MeV/c, 60 < E n < 100.5 MeV and 
12 < R n < 28 cm. 

To increase the statistical power of any observed sig- 
nal events (Section HIE), a tighter kinematic region was 
defined as a subset of the standard region to further sup- 
press K^2 and K e 4 backgrounds. As shown in Figure 14, 
K V 2 and K e 4 events were not uniformly distributed in 
the signal region. The K^i target-scatter events were 
more uniformly distributed in the signal region except 
near the K„2 peak. The imposition of the t:vD{\+2) trig- 
ger on the K e 4 momentum distribution shown in Figure 1 
caused the K e4 background to peak around 160 MeV/c as 
described in Section IIIC2. The accepted K + —> n + vv 
spectrum was monotonically increasing with momentum 
in the signal region (Figure 17). Based on these obser- 
vations, the kinematic region that maximized the ratio 
of the signal acceptance to the total K^i and K e 4 back- 
ground was 165 < P T < 197 MeV/c, 72 < E n < 100 MeV 
and 17 < R n < 28 cm. 

We also defined the "if^-peak region" by the require- 
ments 199 < P n < 215 MeV/c, 100.5 < E n < 115 MeV 
and 28 < R n < 35 cm (Figure 4). Events in the -fC^-peak 
region were employed to set selection criteria, estimate 
background and determine the signal acceptance. 

5. Muon identification 

Muon backgrounds were rejected largely based upon 
the positive identification of the ir + in the range stack 
by the observation of the it — > fi — > e decay chain and by 
the range-momentum relationship. 

The 7T + identification algorithms of -kvv(1) [2] were 
adopted for this analysis, and only a brief description is 
provided here. The analysis of the waveform provided 
by the transient digitizers (TDs) was used to identify the 
7r + — > n + v^ decay in the range stack element that con- 
tained the stopping pion. A neural network was trained 
using kinematically identified 7r + and /i + that stopped in 
the range stack [2]. The [i + — ► e + v e D fl decay was iden- 
tified by TDC information in the range stack counters 
near the stopping counter. 



Since the previous Ttvi>(2) analyses [35, 45] had shown 
that muon backgrounds were small, less restrictive re- 
quirements on the 7r — > /i — > e decay than those in the 
m>D(l) analysis were used for the standard twv{2) re- 
quirements. This provided a 10% increase in the sig- 
nal acceptance. A muon rejection of 133.0 ± 10.7 (Sec- 
tion III C 3) was obtained with a looser requirement on 
the neural network output and no identification of the 
ji — > e decay. The m/&(l) requirements were used to de- 
fine a tighter cut that was used to subdivide the signal 
region as described in Section III E. The tighter cuts had 
a muon rejection of 409.1±60.9. In the following we refer 
to the 7r — > (i — > e requirements as the "TD" cuts. 

The ability to separate pions from muons using the 
range and momentum measurements can be seen in Fig- 
ure 4. The separation was based on the "RNGMOM" cut 
that was placed on the quantity Xrm = {R rs — Rutc) / 
where Rjjtc {vr) was the expected range (uncertainty 
in range) for a given ir + momentum and R rs was the 
measured range in the range stack. 



6. Delayed coincidence requirements 

Determining that the incoming K + came to rest in 
the target was accomplished by observing the delay be- 
tween the incoming particle and the outgoing charged 
track. This requirement rejected incoming beam pions 
that scattered in the target as well as the products of 
K + decay-in-flight. The delayed coincidence also served 
to suppress the CEX background. 

For the standard delayed coincidence requirement, the 
average time of the kaon fiber hits (tx) had to be at least 
3 ns earlier than the average time of the pion fiber hits 
(tjr). The previous nuu(2) analyses [35, 45] used a tighter 
requirement oit^—tx > 6 ns. The looser requirement 
in this analysis resulted in a 9% relative acceptance in- 
crease. Since under certain conditions the resolution on 
tx or t n was degraded, the degraded time resolution was 
taken into account by increasing the minimum delayed 
coincidence allowed. It was increased to 4 ns when the 
energy deposit in the target kaon fibers was less than 50 
MeV, when the time of any of the individual kaon fiber 
differed from the average kaon fiber time by more than 
2 ns or when the time of an individual target pion fiber 
differed from the average time of the pion fiber hits by 
more than 3.5 ns. It was increased to 5 ns when the dif- 
ference between tx and the B4 hodoscope hit time was 
greater than 1 ns or when t n was determined using IC 
hit time because the average time of the pion fiber hits 
had poor resolution. It was increased to 6 ns when the 
difference between t n and t rs was greater than 1.5 ns. 

A tight version of the delayed coincidence with the 
requirement of t n — tx > 6 ns was used to subdivide the 
signal region as described in Section HIE. 
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7. Photon veto requirements 

An event was rejected by the photon veto cut when 
the total energy in a sub-detector within a time window 
exceeded a given threshold. The time window was ref- 
erenced to t rs , the reconstructed time of the pion in the 
range stack. The time window and energy threshold was 
set for each sub-detector using an optimization algorithm 
described in [2] . The end caps were treated as three sep- 
arate sub-detectors EC ou t er , ECi nner and EC 2n d in the 
optimization. ECi nner was the inner ring of the upstream 
EC and had higher accidental rates than the remaining 
EC elements due to its proximity to the incoming beam. 
EC-ind was the EC energy identified by a double-pulse- 
finding algorithm using CCD information. EC ou ter com- 
prised the EC elements not included in ECi nner . The 
optimization procedure determined the rejection and ac- 
ceptance as the time window and energy threshold were 
varied. The optimization goal was to maximize rejection 
for a given value of acceptance. The acceptance sam- 
ple used by the optimization procedure was derived from 
K^2 monitor trigger data. 

The photon veto requirements for the ttuD{1) analysis 
were optimized using K^i peak events that were the dom- 
inant background with photons. Ideally the ■kuv{2) pho- 
ton veto requirements would have been optimized on a 
sample of K W 2 target-scatter events; however, given that 
photon veto rejection needed to be 0(2500), we were un- 
able to prepare such a sample with sufficient statistics, 
0(250000) events, needed to minimize bias in the opti- 
mization result. In lieu of this sample, we optimized the 
photon veto requirements for a majority of sub-detectors 
using a sample of K^i peak events and then optimized the 
requirements for the remaining sub-detectors using mul- 
tiple samples of K n2 target-scatter events as described 
below. 

The main sample of K n 2 target-scatter events failed 
either the CCDPUL cut (Section IIIB2) or the Beam 
Likelihood cut (Section IIIB3) and contained 26317 and 
52621 events in the 1/3 and 2/3 data samples, respec- 
tively. Other K„2 target-scatter samples were composed 
of events failing these cuts or the other target cuts de- 
scribed in Section IIIB2. The size of the other samples 
ranged from 11037 (22037) to 29899 (59871) in the 1/3 
(2/3) data samples. These samples overlapped one an- 
other, but they contained pions with different relative 
populations of the pion scattering angle with respect to 
the beam direction. 

An additional sample, dubbed the "kink" sample con- 
taining K V 2 target scatters where the tt + track had an 
identifiable kink in the xy projection, was created by pro- 
cessing every ■ni>v{l+2) event with a kinked-track recon- 
struction algorithm. For kink reconstruction, the restric- 
tions on the pion fiber energy were removed as well as 
the requirement that the pion fibers had to be within 1 
cm of the extrapolated UTC track. The following criteria 
defined a valid kink event: (1) the event had at least two 
pion fibers that deviated from the UTC extrapolation, 



TABLE I: Time window and energy threshold of the primary 
and secondary photon veto requirements for each sub-detector 
as described in the text. The time window was defined with 
respect to t rs , the reconstructed time of the n + in the range 
stack. RS and TG label the range stack and target parame- 
ters, respectively. The parameters for the sub-detectors below 
the double line were optimized separately as described in the 
text. BVLsame had the additional requirement that recon- 
structed z position satisfy \z\ < 4 cm. 





Primary 




Secondary 




Time Threshold 


Time 


Threshold 


Sub-detector window (ns) 


(MeV) 


window (ns) 


(MeV) 


BV 


[-5.7,10.2] 


0.2 


[-7.5,10.2] 


0.7 


BVL 


[-4.4,10.7] 


0.3 


[-3.5,10.6] 


0.3 


RS 


[-4.3,4.4] 


0.3 


[-3.3,7.8] 


0.6 


ECouter 


[-4.4,8.0] 


0.4 


[-6.0,9.5] 


0.2 


ECinner 


[-3.7,5.6] 


0.2 


[-14.0,9.1] 


0.2 


EC'2nd 


[-5.7,2.5] 


10.6 


[-5.7,2.7] 


0.2 


TG 


[-2.7,2.2] 


2.0 


[-6.6,2.3] 


1.7 


IC 


[-2.0,4.5] 


5.0 


[-2.9,9.3] 


1.4 


vc 


[-6.6,1.8] 


6.8 


[-7.5,7.1] 


5.0 


CO 


[-0.1,5.9] 


0.6 


[-0.8,5.1] 


6.0 


nco 


[-5.5,2.3] 


3.0 


[-4.5,3.3] 


1.6 




[-35.7,-5.7] 


30.0 


[-37.5,-7.5] 


30.0 


AD 


[-2.0,8.0] 


0.6 


[-2.0,8.0] 


0.6 


DPV 


[-5.0,10.0] 


0.0 


[-5.0,10.0] 


0.0 


BV L same 


[-5.0,-2.0] 


10.0 


[-5.0,-2.0] 


10.0 



(2) at least one of the fibers from (1) must be adjacent 
to a kaon fiber, (3) the remaining pion fibers must be 
along the UTC extrapolation and (4) the event must be 
rejected by the criteria placed on the standard target- 
track reconstruction. The final criterion guaranteed that 
the kink sample was independent of the sample of signal 
events and the other samples described in the previous 
paragraph. Although the resulting kink sample had only 
11833 events, it provided a sample rich in target scatters 
that was used in understanding the response of the AD 
as described below. 

Before beginning the photon veto optimization proce- 
dure, we applied a cut on the activity in the BV prior 
to t rs (BV ear i y ) because a large energy deposit (> 30 
MeV) in the BV prior to the kaon decay would pre- 
vent the TDCs from registering activity coincident with 
t rs [47]. The -kvv(1) set of parameters as listed in Ta- 
ble VI of [2] was the starting point for the nvv[2) opti- 
mization that included all sub-detectors except the AD 
and DPV. Primary and secondary sets of parameters 
were determined for the eleven sub-detectors listed in Ta- 
ble I. 

The AD parameters were determined using the kink 
sample after application of a photon veto cut with a re- 
jection of approximately 175 with looser settings on the 
parameters of the other sub-detectors. These AD pa- 
rameters yielded an additional photon veto rejection on 
the main K„2 target-scatter sample of 1.95 ± 0.08 with a 
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8. Single- cut failure study 
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FIG. 10: The offline rejection vs. total acceptance for the 
optimized photon veto cuts. The error bars represent the sta- 
tistical uncertainty. The labeled starred points are described 
in the text. 



94% acceptance factor [47] . The main target-scatter 
sample was also used to optimize the DPV parameters. 
After application of all other photon veto parameters at 
the primary setting listed in Table I, the DPV rejection 
was measured to be 1.13 ±0.09 with an acceptance factor 
of 99.99%. 



The BVL same cut was devised subsequent to the 
single-cut failure study on the 1/3 sample (Sec- 
tion III B 8) . The cut removed potential K V 2 background 
when both photons from the 7r° decay deposited energy 
in the same BVL element. Such an occurrence yielded 
a reconstructed time earlier than t rs , a reconstructed z 
position near the center of the element and an apparent 
energy greater than 10 MeV. 



Figure 10 shows the offline rejection for fixed values of 
of the total (online and offline) acceptance for the photon 
veto. The parameters in the primary column in Table I 
corresponded to the standard photon veto cut ("Stan- 
dard" in the Figure) . For the more restrictive ( "Tight" ) 
photon veto cut described in Section HIE, events were 
rejected that failed the criteria established by either the 
primary or the secondary parameters. The additional 
settings labeled "Loose" and 'WP(l)" in the Figure, of 
the photon veto cuts were used for background estimation 
(Section III C 3) and consistency checks (Section III C 9). 



After determination of the selection criteria using the 
1/3 data sample, we performed a "single-cut" failure 
study to identify unexpected sources of background or 
potential analysis flaws. Individual cuts that exploited 
similar background characteristics were grouped together 
to form the following twelve cut categories: 

1. The cuts on R w , P n and E v (Section HIB4). 

2. All photon veto cuts except those on the AD and 
target. 

3. The photon veto cut on the AD. 

4. The target photon veto and OPSVETO (Sec- 
tion IIIB2) cut. 

5. The delayed coincidence cut (Section IIIB6). 

6. The ir/n range-momentum separation requirement 
(Section III C 3) and the pion track requirements 
(Section III B 3) excluding the cuts in the next two 
categories. 

7. The Beam Likelihood cut (Section IIIB3). 

8. The cut requiring the z position of the extrapo- 
lated UTC track to be more than 6.5 cm from the 
upstream end of the target (Section IIIB3). 

9. The CCDPUL and kaon fiber timing cuts (Sec- 
tion IIIB2). 

10. The cuts related to the identification of the tt — > 
[i — ► e decay chain. 

11. The cuts to suppress beam-related backgrounds. 

12. The cuts on pion fiber energy, the pattern of kaon 
and pion fibers and the target-track fitter (Sec- 
tion IIIB2). 

All events in the 1/3 data sample that failed only one of 
these twelve categories were examined. We found four 
events that contained evidence of two potential analysis 
flaws. 

Three of the events not rejected by the photon veto 
cuts showed evidence of a large energy deposit in the 
BVL. These events were shown to be due to iV-2 decays 
in which both photons from the 7r° decay deposited en- 
ergy in the same BVL counter [48]. The simultaneous 
activity at each end of a BVL element led to an erro- 
neous average time prior to t rs that was outside the veto 
time window. The u BVL same v cut, previously described 
in Section III B 7, was devised to remove these events. 

The remaining event of the four failed only the pho- 
ton veto criteria in the AD and revealed a potential flaw 
in the CCDPUL target-pulse fitting algorithm when the 
fitted time of the first pulse was inconsistent with the 
average kaon fiber time. The inconsistency arose when 
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FIG. 11: Schematic representation of a K^z target-scatter in 
which the 7r + initially traveled in the z direction, scattered in 
a kaon fiber and was directed into the barrel region. The two 
photons from the decay of the recoiling it were directed into 
the beam region. 

the fitting algorithm incorrectly assigned the first pulse 
time to an actual second pulse because the second pulse 
energy was larger than the first pulse energy. The CCD- 
PUL timing criteria already described in Section III B 2 
were developed to remove the analysis flaw. 

No analysis flaws or unexpected sources of background 
were revealed by the "single-cut" failure study of the 2/3 
data sample. 



C. Backgrounds 

1. K n2 -related background 

The if^-related background contained three compo- 
nents: Kn2 target-scatter, K^i range-stack-scatter and 
Kn2y. In order to have a K W 2 event in the signal region, 
the photons from the ir° decay had to escape detection 
and the 7r + had to lose energy via scattering or resolution 
effects such that it fell into the signal phase space. Scat- 
tering could happen in the target [K^i target-scatter) or 
in the range stack (K^ range-stack-scatter). The target- 
scatter component dominated in the background. Since 
the 7r + from a K n %^ decay was not monochromatic, the 
three final-state photons needed to escape detection for 
this type of event to be a background. 

The topology of the most problematic type of 
target-scatter was that of a 7r + initially traveling in the 
kaon fibers and scattering into the barrel region of the 
detector as indicated schematically in Figure 11. This 
type of target-scatter was difficult to reject because some 
energy deposited in the target by the scattering ir + oc- 
curred in a kaon fiber (Section III B 2) and could not 
always be distinguished from the larger energy deposited 
by the stopping kaon (Figure 7). In addition the 7r° was 
also traveling parallel to the beam direction and the re- 
sulting photons from the 7r° decay were directed at the 
upstream or downstream ends of the detector where the 



photon veto was less efficient. 

In the K^2 target-scatter background estimate, the 
two bifurcation cuts chosen were: the standard photon 
veto cuts (CUT1) and the target-quality cuts (CUT2), 
since both of these gave powerful rejection of the 
target-scatter background. The bifurcation analysis sam- 
ple was prepared by applying cuts to remove the contam- 
ination from muon, beam and charge- exchange events. In 
the normalization branch, a combination of K V 2 target- 
scatter, K^2 range-stack-scatter and K^-y events were se- 
lected by inverting the photon veto cuts (CUT1). All the 
target-quality cuts (CUT2) were applied to the sample, 
resulting in 1131 events left in the normalization branch 
in the 2/3 sample. After corrections for range- 
stack-scatter contamination (detailed below), A^ norm = 
1107.7±33.8(stat.)l2's( s y s t-) events remained in the nor- 
malization branch. The systematic uncertainty is due to 
the correction for range-stack-scatters. Corrections for 
contamination due to K^y are discussed later in this 
Section. 

For the K r 2 target-scatter rejection branch, the 
K^2 target-scatter events were classified into two non- 
exclusive categories. The first category, "z-scatter" , oc- 
curred when the 7r + traveling parallel or anti-parallel to 
the beam direction scattered in a kaon fiber into the 
barrel region of the detector as depicted in Figure 11. 
The second category, "xy-scatter" , occurred when the 
7T + scattered outside of the kaon fibers, and the scatter 
was visible in the xy plane. To measure the rejection 
of the photon veto for target-scatter events, six classes 
of K^2 target-scatter events containing varying mixtures 
of xy-scatter and z-scatter events were created by ap- 
plying or inverting various combinations of the require- 
ments on 7r + in the target (Section IIIB 2). The primary 
K^2 target-scatter sample (Section III B 7) , considered 
to be the richest in z-scatters, was chosen to measure 
the photon veto rejection, giving 52621 events for the re- 
gion C + D (Figure 5). The photon veto cuts (CUT1) 
were then applied to the remaining K„2 target-scatter 
events, leaving 22 events for the region C for a rejection 
of 52621/22 = 2392±510, where the uncertainty is statis- 
tical only. The pion momentum distributions of the nor- 
malization and rejection branches are shown in Figure 12. 
The photon veto rejections measured in the other five 
classes of target-scatter events were consistent with that 
of the primary sample and extended from a minimum of 
2204 ±697 to a maximum of 2758 ±650. The photon veto 
rejection of the main K^2 target-scatter sample in the 
tighter kinematic region (Section III B 4) was measured 
to be 2193 ± 517 in agreement with the measurement in 
the standard kinematic region. The range of measured 
rejection values in the five other classes was used to set 
the systematic uncertainty in the photon veto rejection 
on the K V 2 target-scatters, giving a photon veto rejec- 
tion of 2392 ± 510±f||. Using Equation (3), the uncor- 
rected number of K W 2 target-scatter background events 
was measured to be 

b% = 3/2 x (1107.7±33.8j£j)/((2392±510t?|f) - 1) 
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FIG. 12: (Top) Pion momentum distribution for the normal- 
ization branch in the 2/3 data sample. The black histogram 
represents the distribution after inversion of the photon veto 
cut. The blue dashed histogram shows the suppression of 
all the target cuts except for the CCDPUL cut and the red 
shaded histogram shows the suppression of all the target cuts 
including the CCDPUL cut. (Bottom) Pion momentum dis- 
tribution of the rejection branch in the 2/3 data sample. 
The events obtained from the inversion of the target cuts are 
shown in black and the red shaded histogram shows the events 
that survive the photon veto cut. The events surviving the 
photon veto cut that peak near 236 MeV/c were due to K^2 
decays. The cuts on pion kinetic energy and range were not 
applied for these distributions in order to show the momen- 
tum distribution outside the ttuv{2) signal region. 



= 0.695 ±0.150™>, (4) 

where the first uncertainty was statistical and the second 
uncertainty systematic. 

For the range-stack-scatter background events, 

the cuts with the most powerful rejection were the range- 
stack track quality and the photon veto cuts. The K V 2 
range-stack-scatter normalization branch was a modified 
version of the K^i target-scatter normalization branch, 
with the range-stack quality cuts inverted instead of be- 
ing applied before the inversion of the photon veto cut as 
was done in the target-scatter normalization branch. 
This sample of N 2 = 281 events was heavily contami- 
nated with target-scatter events due to the inefficiency 
of the range-stack-scatter cuts. The Ni = 1131 events 
remaining at the end of the K V 2 target-scatter normal- 
ization branch consisted of N tg target-scatter events with 
contamination due to N rs range-stack-scatter events: 

N ta +N rs =Ni. (5) 

These N tg K n2 target-scatter events and N rs K V 2 range- 
stack-scatter events were also related to the N2 events 
remaining in the range-stack-scatter normalization 
branch by 

x N tg + (R rs — 1) x N rs = N 2 , (6) 

where A rs = 0.888 ± 0.001(stat.) ± 0.012(syst.) was the 
acceptance factor for the range-stack quality cuts and 
R rs = 7.06 ± 0.47 was the rejection of K V 2 range-stack- 
scatter events by the range-stack quality cuts measured 
using events with momentum consistent with the K v2 - 
peak region, but range and energy in the irvD(2) signal 
region as would be expected for a range-stack-scatter. 
The systematic uncertainty on A rs was due to the larger 
uncertainty on the measured kinematics of the 7r sca t mon- 
itor data used to measure the acceptance of the range- 
stack quality cuts as described in Section HID. By solv- 
ing Equations (5) and (6) simultaneously, it was possi- 
ble to estimate the number of K n2 target-scatter events 
(N t9 = 1107.7 ±33.8±||) and the number of K n2 range- 
stack-scatter events (N rs — 23.3±3.5lg'Q) present in the 
original K n2 target-scatter normalization branch, where 
the first uncertainty is statistical and the second is sys- 
tematic due to the acceptance factor A rs . 

The photon veto rejection on the K n2 range-stack- 
scatter events should be the same as that for the unseat - 
tered peak events as the back-to-back correlation of 
the 7r + and ir° was maintained. The K w2 range-stack- 
scatter rejection branch was created by applying all cuts 
other than the photon veto cuts and the pion kinematic 
cuts (Section III B 4) . The K„ 2 peak region events were 
selected, creating a sample of 122581 events for the region 
C + D. The PV cuts (CUT1) were then applied to the 
remaining K v2 events, leaving 106 events in region C for 
a photon veto rejection of 122581/106 = 1156±112. The 
number of K^2 range-stack-scatter background events 
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was measured to be 



3/2 x (23.3 ± 3.5±|;g)/((1156 ± 112) 
0.030 ± 0.005(stat.) ± 0.004(syst.). 



1) 



(7) 



The Kn2~/ background estimate did not use the bifur- 
cation method, but used a combination of K v2 events se- 
lected in the irvv{l+2) data and simulated K V 2 and K^j 
events. We simulated both the inner brcmsstrahlung 
(dominant) and direct emission amplitudes of K^-y de- 
cays assuming no interference between them [49]. The 
K V 2-y background was estimated as 



N 



(8) 



where 



N = 106 was the number of isT^-peak events in the 
nvD (1+2) trigger sample after all selection criteria 
were applied except that the if^-peak region kine- 
matic region was selected. 

k = 417 ± 24 was the ratio of acceptance factors of 
Kn2 to -fC^ events as determined from simulated 
K n 2 and K n 2-y decays taking into account the 
and K W 2-y branching ratios. We used B(K + — > 
7r + 7r°7) = (11. 1±0. 6) x 10~ 4 obtained by correcting 
the partial branching ratio (2.75 ±0.15) x 10~ 4 [24] 
measured for 55 < < 90 MeV to the full energy 
range using simulation. 



i? 7 = 5.04 ± 0.10 was the additional photon veto re- 
jection afforded by the radiative photon. This ad- 
ditional rejection was calculated by combining the 
distribution of the radiative photon from simulated 
events with the measured single photon detection 
efficiency as a function of angle and energy from 
K n2 data [48]. 

The final anticipated number of K^i background events 
was V 1 = 0.076±0.007±0.006, where the first uncertainty 
was statistical and the second was systematic (due to n 
and Fly). 

The inverted photon veto used to select events for the 
K V 2 target-scatter normalization branch would have also 
selected K n 2-y events. We corrected the estimate for the 
K V 2 target-scatter background by subtracting 6 7 , 

fjtg = tfg _ b i 

= 0.619 ±0.150tHoo • 



2. 



K+ — > 7r + 7r e + v background 



Despite the small branching ratio of (4.09 ± 0.10) x 
10~ 5 [24], K e4: could be a background if the tt~ and the 
e + escaped detection in the target. The distribution of 
the sum of the kinetic energies of the 7r _ and the e + 
vs. the reconstructed 7r + momentum in simulated events 
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FIG. 13: The sum of the n~ and e kinetic energies versus the 
reconstructed momentum of the ir + for simulated K e 4 events 
that passed the ttvv(1+2) trigger. 



shown in Figure 13 indicates where the K ei background 
would occur kinematically. 

Since the main characteristic of K e 4 event was extra 
energy in the target from the n~ and the e + , the tar- 
get photon veto (TGPV), OPSVETO and CCDPUL cuts 
were the most effective cuts to suppress this background. 
Due to contamination by other types of background, such 
as if W 2 _ target-scatter, it was not possible to isolate a pure 
K e 4 background sample for a bifurcation analysis using 
data only. Nonetheless, a if e4 -rich sample was selected 
from data using the CCDPUL • TGPV • OPSVETO re- 
quirement and served as the normalization branch. We 
established that the majority of the events in the nor- 
malization branch were likely to be due to K e ± decays by 
removing the CCDPUL requirement and comparing the 
momentum distribution of the selected events in the 1/3 
sample with the expectation from simulation (Figure 14). 
In addition, the target information for the 69 events in 
Figure 14 was visually examined and the events clas- 
sified based on the topology, ionization pattern, curva- 
ture, range and energy of the putative tracks. Fifty-nine 
events were classified as K e 4, three as K + — » n + fi~(i + 
(9) and four as K v2 or K^ 2 -f (including events with apparent 
7T° — > e + e~7 decays). The classification of the remaining 
three events was ambiguous. Assuming half the ambigu- 
ous events were K e 4 yielded a purity of <~ 88%. One 
example of a K e 4 candidate event is shown in Figure 6. 

Simulated K e 4 events were used to determine the re- 
jection of the TGPV, OPSVETO, and CCDPUL require- 
ments. Negative pion absorption in the target was mod- 
eled based on the energy spectrum of stopped 7r~ in 
plastic scintillator observed in E787 [50]. We assumed 
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FIG. 14: Momentum distribution of the tt + in the K e 4 nor- 
malization branch for the 1/3 sample before the application 
of the CCDPUL cut. The points represent the data. The 
unshaded histogram is the distribution as expected from sim- 
ulated K e 4 events. The histogram area is normalized to the 
number of data events. The shaded histogram represents the 
normalization branch for K„z target-scatter events. The ratio 
of the area of the shaded to the unshaded histogram has been 
arbitrarily set to 1/8 times the ratio of the target-scatter 
background to the K e 4 background for display purposes. 



that all energy generated from ir~ absorption would be 
promptly deposited in the single fiber where the n~ came 
to rest. This assumption conservatively neglected the 
possibility that detectable activity from tt~ absorption 
could occur elsewhere in the detector 1 . Positron in- 
teractions were well-modeled in our EGS4-based simu- 
lation [46]. The rejection of the CCDPUL, TGPV and 
OPSVETO requirements were correlated because the tar- 
get fibers containing the deposited energy of the tt~ and 
e + could have been classified as kaon, pion, photon or 
opposite-side pion fibers. We used the energy of the sim- 



Negative pions are predominantly absorbed by carbon nuclei 
in scintillator resulting in multi-nucleon emission. The mea- 
sured [51, 52] rates of emission per stopped n~ of neutrons, 
protons, deuterons, tritons and alphas are approximately 2.8, 
0.3, 0.2, 0.1 and 0.6, respectively, with typical kinetic energies 
of tens of MeV. Because of their short range, these charged par- 
ticles will deposit energy very close to their absorption points. 
The mean free path of emitted neutrons is tens of cm leading 
to energy deposition relatively far from the absorption point. In 
addition, the residual nucleus is unstable and can deexcite by 
emission of photons with typical energies of 1-2 MeV. The en- 
ergy spectrum measurement [50] was sensitive to energy near 
the absorption point and thus largely neglected any additional 
energy deposition due to the latter two processes. 



ulated deposits to estimate the rejection of these cuts as 
R = 52+2,g 1 . As we did not precisely model either n~ ab- 
sorption or the inactive material of the target such as the 
gaps between the fibers and the cladding and wrapping 
material of each fiber, we varied the threshold for the 
energy treated by the CCDPUL (TGPV • OPSVETO) 
cut by a factor of 5 (1.5) to estimate the systematic 
uncertainty associated with the rejection of these cuts. 
The normalization branch in the 2/3 sample contained 
6 events so the K e 4 background was measured to be 
3/2 x 6/(521^ - 1) = 0.176 ±0.072±°;^| events, where 
the first error was statistical and the second was system- 
atic. 



3. Muon background 

The decays K + — > fi + u^, K + — > /i + f M 7 and K + — » 
/i + 7T could contribute background in the nvu(2) kine- 
matic region as indicated in Figure 4. All three processes 
required the muon to be mis-identified as a pion in order 
to be a background. The first decay would be background 
if the kinematics of the fi + were mis-reconstructed and 
the latter two decays would be background if the photons 
went undetected. 

The two bifurcation cuts were the tt — ► ji — > e iden- 
tification or "TD" cut (CUT1), and the 7r//i range- 
momentum separation or "RNGMOM" (CUT2). The 
normalization branch defined by inverting the TD cut 
yielded zero events in the 2/3 sample, so iV norm was as- 
signed to be 1 event. 

The fi + rejection branch contained C + D = 20488 
events in the 2/3 sample and was selected by inverting 
CUT2 and applying cuts to remove beam backgrounds 
and the -nvviV) version of the photon veto cut (Figure 10) 
to suppress K n2 backgrounds. After the application of 
the TD cut, the number of events remaining was C = 154 
for a measured TD cut rejection of 133.0±10.7. Thus, the 
background was estimated to be 3/2 x (1±1)/((133.0± 
10.7) - 1) = 0.0114 ±0.0114. 

4- Charge- exchange background 

When the K° due to CEX in the target decayed as a 
K\ it was a potential background. The delayed coinci- 
dence requirement effectively removed any contribution 
from the short-lived Kg. The semileptonic decay pro- 
cesses — > 7r + e~i? e and K Q L — > tt + yTv^ with branching 
ratios of 20% and 14%, respectively, were considered to 
be the most likely to form a background. 

The CEX background could also contain a component 
due to hyperon production where a 7r + was either pro- 
duced with the hyperon or was a hyperon decay product. 
Hyperon production would result from X°-nucleon inter- 
actions if the K° oscillated to a K°. 

Simulation studies showed that there was often a gap 
between the pion and kaon fibers and that the recon- 
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FIG. 15: Schematic diagram of the charge- exchange process 
in the target. 



5. Single-beam background 

The bifurcation cuts for single-beam background were 
the delayed coincidence (CUT1) and B4 energy of less 
than 1.0 MeV (CUT2). CUT2 selected beam pions. The 
sample was selected by applying all the photon veto cuts 
except TGPV, the kinematic cuts, TD cuts and beam 
cuts except the delayed coincidence and CUT2. The 
normalization sample formed by the inversion of the de- 
layed coincidence cut samples yielded zero events, so that 
N nolm was assigned to be 1. The rejection sample con- 
tained C + D = 12850 events in the 2/3 data. After the 
application of the delayed coincidence, C — 2 events re- 
mained for a rejection of the delayed coincidence cut of 
6425 ±4543. The measured single-beam background was 
3/2 x (1 ± 1) /(6425 - 1) = 0.00023 ± 0.00023 events. 



structed z of the pion track was not consistent with the 
energy deposited in the kaon fibers as indicated schemat- 
ically in Figure 15. A CEX-rich sample that served as 
the normalization branch was selected in 7rw(l+2) trig- 
gers by requiring a gap between the pion and kaon fibers. 
No target energy cuts were applied in the selection of the 
normalization sample since the lepton from K® L decay or 
the ir~ and ir° associated with hyperon production might 
deposit extra energy in the target. The offline delayed co- 
incidence requirement was also not applied for the nor- 
malization sample. In contrast to the previous ttv9(2) 
analysis [45], the ability to create a normalization sam- 
ple from ttv9(1+2) data avoided uncertainties associated 
with the effective CEX cross-section and related efficien- 
cies. In addition a normalization sample selected from 
the data contains contributions from all CEX processes 
including and hyperon decays. The normalization 
branch in the 2/3 sample contained one event. 

The rejection associated with finding a gap in the CEX 
events, the target energy cuts and the delayed coinci- 
dence was determined from simulated CEX — > tt + £~u 
events. For the simulation of CEX events, reconstructed 
Kg — * 7r + 7r~ events obtained from the CEX monitor trig- 
ger data (Section II B) were used as the K\ production 
point and momentum. The rejection of the delayed coin- 
cidence and gap-finding cuts exploit the flight of the neu- 
tral kaon and should be the same for K\ — > ir + £~D and 
hyperons. The rejection of the target energy cuts on the 
7r~ or 7T° associated with hyperon decay was estimated to 
be comparable to or greater than the rejection of the tar- 
get energy cuts on the lepton from K® — > ir + £~i> decays. 
The background was measured in the 2/3 sample to be 
0.013 ± 0.013(stat.)j;°;°J°(syst.) events. The systematic 
uncertainty was estimated by the same threshold varia- 
tions of the target energy cuts that were used for K e 4 
(Section III C 2). 



6. Double-beam background 

Double-beam background had two components, KK 
and Kir. For the KK (Kit) background, the decay prod- 
ucts of the initial kaon were undetected and a subsequent 
kaon decay (a scattered beam pion) provided the outgo- 
ing 7T + . 

The KK component of the background b KK was deter- 
mined using the bifurcation procedure described in Sec- 
tion III A 3 with modifications to compensate for poor 
statistics: 

b KK = fx n KK /TKK (10) 

Rkk — 1 

where / = 3/2 was the scale factor for the 2/3 data 
sample (Equation (3)). 

Rkk = 1576/4 = 394 ± 197 was the measured rejection 
of the kaon Cerenkov and BWPC cuts (CUT1) on 
a beam kaon at the time of the outgoing pion (t rs ). 
The rejection sample was prepared by vetoing beam 
pions at t rs via Cerenkov information, by requiring 
a second track at t rs in the B4 counter with an 
energy deposit consistent with a kaon and by the 
more stringent target gap requirement described in 
Section IIIB2. The latter criterion ensured activ- 
ity in the target in two spatially and temporally 
distinct regions indicative of a double-beam event. 

nKK was the normalization provided by a second bifur- 

r KK i . . 

cation of the standard normalization branch. The 
second bifurcation [2] exploited the lack of correla- 
tion between the AD and target cuts to improve the 
statistical power of the measurement. The normal- 
ization branch was prepared by inverting CUT1, by 
vetoing entering pions at t rs using Cerenkov infor- 
mation and by application of the n — > /i — > e and 
track quality requirements on the outgoing track 
and contained 2699 events. The application of the 
AD photon veto cuts reduced the sample to 325 



19 



events for a rejection of Tkk — 8.3 ± 0.4. The 
application of the target cuts on the 2699 events 
yielded zero events so we assigned tlkk = 1- 



lKK 



These values gave the KK background of b 
0.00046 ± 0.00046 events. 

An analogous method was used to estimate the Kir 
component of the double-beam background 



1 - 0.606 



x / x 



R 



Kit 



1 



(11) 



where the additional scale factor of 1/(1 — 0.606) was 
included to correct for the data accumulated with the 
online pion Cerenkov veto in the iru9(2) trigger (Sec- 
tion II B). With the online veto, the offline rejection of 
the pion Cerenkov cuts was low and the normalization 
branch lacked statistics. To obtain b K7r , the Kir back- 
ground estimate obtained without the online veto was 
scaled by the ratio of the kaon exposures. As a check 
of this scaling, we verified that the KK background es- 
timates for the two trigger configurations were consis- 
tent and that the Kir background for the iru9(l) trigger 
data was consistent for the two data-taking periods [53] . 
With the measured values of R K ^ = 2467/4 = 617 ±308, 
n & = 1 ± 1 and r K7T = 4435/464 = 9.6 ± 0.4 in Equa- 
tion (11), we obtained b K7T = 0.00064 ± 0.00064 events. 



7. Background summary 

The contribution of each background component is 
listed in Table II. The total background was estimated 



Process 


Background events 


Ktt2 target-scatter 
K n 2 range-stack-scatter 

K 7T 2~f 

K e4 
Muon 

Charge-exchange 
Beam 


0.619 ±0.150to'ioo 
0.030 ± 0.005 ± 0.004 
0.076 ± 0.007 ± 0.006 
0.176 ± 0.072±g-?H 
0.011 ± 0.011 
0.013 ± 0.013±g"go§ 
0.001 ±0.001 


Total 


0.927 ± 0.168+ 0.237 



TABLE II: Summary of the estimated number of events in 
the signal region from each background component. Each 
component is described in the text. 



to be 0.927± 0.168±g;|l? events and was dominated by 
the K n 2 target-scatter component that was the largest 
contribution to the statistical uncertainty. The system- 
atic uncertainty was dominated by the contribution from 
the -ftT e 4 background due to the inability to establish a 
precise correspondence between the energy observed in 
the target in data and simulation. 

A number of background consistency and validity 
checks were performed as described below. 



8. Background contamination evaluation 

Due to the difficulty of isolating background samples, 
studies were performed to estimate the degree of con- 
tamination (i.e., events due to background from other 
sources) in the K^ target-scatter normalization and re- 
jection branches (Figure 12). 

The effect of muon contamination of the K^i back- 
ground estimate was determined separately for the nor- 
malization and rejection branches with and without the 
ir — > /i — > e (TD) and ir//j, range-momentum separation 
(RNGMOM) cuts. The normalization branch used in the 
if 7r2 -scatter study was assumed to be the sum of ir + and 
H + components such that 



A n 



1131 = N* 



norm 



(12) 



When the TD and RNGMOM cuts were not applied, the 
observed number of events in the normalization branch 



was 



n = 12980 



(13) 



where A v (R^) is the acceptance (rejection) factor for 
the combination of the TD and RNGMOM cuts for pi- 
ons (muons). The rejection of the TD cut was evaluated 
as 133.0 ± 10.7 as part of the muon background estimate 
(Section 1110 3). The RNGMOM rejection of 28.3 ± 1.1 
was evaluated using the muon normalization branch for 
a total muon rejection of R^ = 3764 ± 333. The ac- 
ceptance factor A n of the combination of the TD and 
RNGMOM cuts was determined on samples of AT^-peak 
region events that failed different combinations of the 
target-scatter cuts used to assess the uncertainty in the 
photon veto rejection as described in Section III B 7. The 
acceptance factors for these samples, both before and af- 
ter the application of the standard photon veto cut, were 
consistent and yielded A v = 0.809 ± 0.030. These gave 
the muon contamination of the normalization sample of 
A£ orm /A norm = (2.7 ±0.3) x 10~ 3 . Analogous method- 
ology was used to assess the effect of muon contamina- 
tion in the rejection branch. The calculated photon veto 
rejection after correction for muon contamination was 
R n = 2410 ± 518 to be compared with R = 2392 ± 510 
(Section III CI). Inserting these results into the back- 
ground estimate using the bifurcation method (Equa- 
tion (3)) implied that the muon contamination increased 
the Kk2 background (b) estimate by 



b/b* = f 



A norm / N* 
R-ll 1 R-* 
1.010 ±0.002 



which was considered negligible with respect to the esti- 
mated systematic uncertainty. 

A similar treatment limited the overestimate of the 
Kk2 background due to double-beam contamination to 
be < 0.1%. 

The rejection of K e 4 by the photon veto should be 
less than the photon veto rejection of K^ and K^-y 
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in that there were no photons in the final state. Con- 
tamination of the K V 2 rejection sample by K e ^ events 
would therefore reduce the measured photon veto rejec- 
tion. The measured rejection of 39481/18 = 2193 ± 517 
in the tighter kinematic region, that was defined to sup- 
press K e 4 (Section IIIB4), was consistent with the over- 
all rejection of 2392 ± 510 (Section III CI), indicating 
no significant contamination of the K V 2 rejection sam- 
ple by K e 4. The K^2 normalization branch was defined 
by the inversion of the photon veto cuts and the ap- 
plication of the target-quality cuts, including CCDPUL 
and OPSVETO. The K e ^ normalization was prepared 
by app lication of CCDPUL • TGPV • OPSVETO. Since 
TGPV • OPSVETO • CCDPUL was a subset of the 
normalization branch, the contamination of the K n 2 nor- 
malization branch by K e 4 was less than the six events se- 
lected in the K e 4 normalization branch (Section III C 2) 
and hence negligible compared to the 1131 events in the 
K n2 normalization branch (Equation (5)). 



9. Background consistency checks 

The consistency of the background estimate was 
checked in three distinct data regions just outside the 
signal region that were created by loosening the photon 
veto cut and the CCDPUL cut. The region CCDi was 
immediately adjacent to the signal region and contained 
events with a CCDPUL second-pulse energy above the 
standard threshold of 1.25 MeV and below 2.5 MeV. The 
region PV\ was immediately adjacent to the signal region 
and defined by events rejected by the standard photon 
veto and accepted by the loose photon veto cuts (Fig- 
ure 10). The region PV2 was adjacent to PV\ and defined 
by events rejected by the loose photon veto and accepted 
by the ttuu(1) photon veto cuts (Figure 10). The num- 
bers of expected background events in these regions were 
calculated in the same manner as for the signal region. 

Table III shows the number of expected and observed 
events in the three regions as well as the probability of 
the observed number of events or fewer given the ex- 
pectation. The combined probability of 5% for the two 
regions nearest the signal region may have indicated that 
the background was overestimated, but the re-evaluation 
of this combined probability at the lower limit of the sys- 
tematic uncertainties [54] gave 14% for the two closest 
regions which demonstrated that the assigned systematic 
uncertainties were reasonable. 

The assignment of N DOlm = 1 when no events were ob- 
served in the normalization branch (Section III A 3) was 
only made for the muon, single-beam and KK double- 
beam backgrounds. Thus, this assignment could have 
overestimated the total background by, at most, 0.012 
events or 1.3%. 



Region 


N E 


No 


T(N ;N E ) 


Combined 


CCD! 


0.79^5 





0.45 [0.29,0.62] 




PVi 


9.0911;! 


3 


0.02 [0.01,0.05] 


0.05 [0.02,0.14] 


PV 2 


32.4lg 2 1 3 


34 


0.61 [0.05,0.98] 


0.14 [0.01,0.40] 



TABLE III: Comparison of the expected (Ne) and observed 
(No) number of background events in the three regions 
CCDi, PVi, and PV2 outside the signal region. The cen- 
tral value of Ne is given along with the combined statistical 
and systematic uncertainties. P(No; Ne) is the probability of 
observing No or fewer events when Ne events are expected. 
The rightmost column "Combined" gives the probability of 
the combined observation in that region and the region (s) of 
the preceding row(s). The numbers in square brackets are 
the probabilities reevaluated at the upper and lower bounds 
of the uncertainty on Ne [54]. 

D. Acceptance and sensitivity 

We assessed the overall acceptance of all selection cri- 
teria by dividing the criteria into components that could 
be measured separately using monitor triggers or sim- 
ulated data. Simulated data were used to estimate the 
acceptance of the trigger and decay phase space as well as 
to assess the impact of nuclear interactions. The overall 
acceptance was the product of the acceptance factors for 
each component. Correlated cuts were grouped together 
for evaluation. 



1. Acceptance factors from K ^2 events 

K^2 monitor triggers were used to assess the compo- 
nents of the acceptance regarding the kaon beam, the 
charged track, the event topology and the standard pho- 
ton veto. The acceptance factors are listed in Table IV 
and described below. 

TABLE IV: Acceptance factors of the K + — * -k + vv selection 
criteria measured with A" M 2 monitor trigger data. Only sta- 
tistical uncertainties are shown. The product is Ak^ 2 ■ 
Cut Acceptance factor 

Range stack track reconstruction 0.99993 ± 0.00001 
UTC-range stack track matching 0.99943 ± 0.00002 
Beam and target pattern 0.15081 ± 0.00018 

Photon veto 0.48122 ± 0.00200 

A K „ 2 0.07253 ± 0.00031 

To measure the acceptance of the range stack track 
reconstruction, a sample of monitor triggers was se- 
lected by requiring a good track in the target and UTC, 
an energy deposit in the B4 hodoscope consistent with a 
beam kaon and a delayed-coincidence of > 5 ns based on 
Ck and the IC. 

The acceptance associated with the matching of the 
range stack and UTC track was assessed using a 
sample with a good track in the range stack, a delayed- 
coincidence of > 5 ns based on the Ck and the IC, and 
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a single entering kaon selected based on the B4 energy 
deposit and the beam Cerenkov and wire chambers. 

The acceptance factor associated with the beam and 
target pattern recognition was evaluated on a sample 
of K„ 2 events that were required to have a single en- 
tering kaon and a good track in the UTC and range 
stack with |cos#| < 0.5. In addition the momentum 
of the reconstructed track was required to be within two 
standard deviations of the expectation for K^ 2 decays. 
There were over forty individual cuts associated with 
the beam and target pattern recognition as described 
in Section IIIB2. The majority of the individual cuts 
had acceptance greater than 90% except for the delayed- 
coincidence (75.5%) and the CCDPUL (45.1%) require- 
ments. 

To measure the acceptance of the standard photon 
veto, an additional criterion was applied to the K^ 2 
events used for the beam and target acceptance. As 
muons from K^ 2 decay can penetrate into the barrel veto 
liner, the reconstructed track was required to stop before 
the outermost layer of the range stack. The acceptance 
factor given in Table IV evaluated in this manner yielded 
the overall acceptance of both the online and offline pho- 
ton veto cuts as the K^ 2 monitor trigger did not include 
the photon veto. 



2. Acceptance factors from K W 2 events 

The K v2 monitor data were used to assess the accep- 
tance factors associated with charged track reconstruc- 
tion in the UTC and in the target. The acceptance of the 
veto of an additional track in the target (OPSVETO, Sec- 
tion III B 2) was also measured with K n2 monitors. The 
factors are listed in Table V and described below. 

TABLE V: Acceptance factors measured with K W 2 monitor 
trigger data. Uncertainties are statistical only. The product 
of all factors is Ak^ 2 ■ 

Cut Acceptance factor 
UTC reconstruction 0.94345 ± 0.00019 
OPSVETO 0.97417 ± 0.00063 
Target track reconstruction 0.71851 ± 0.00181 
A g , 2 0.6604 ±0.0018 

To measure the acceptance of the UTC reconstruction, 
events from the K n2 monitor trigger were required to 
have a well-reconstructed track in the range stack and 
agreement between the online and offline determination 
of the range stack stopping counter. 

For the measurement of the acceptance of the 
OPSVETO cut, in addition to the requirements de- 
scribed above, the charged track was required to be well- 
reconstructed in the UTC and range stack, identified as 
a 7T+ from K v2 decay based on the measured range and 
momentum as well as the ir — > \x — > e signature in the 
stopping counter and kinematically consistent with the 



pion from a K n2 decay. Cuts were also applied to ensure 
a single kaon entered the target. 

In addition to the requirements described above, the 
OPSVETO and target photon veto cuts were applied to 
the K v2 monitor events to assess the cumulative accep- 
tance of the ten cuts associated with target track recon- 
struction. These ten cuts were designed to reject tracks 
that contained an indication of a kink or discontinuity 
in the pattern of target fibers or target fibers with an 
unexpected energy deposit (Section IIIB2). Two indi- 
vidual cuts with less than 90% acceptance were the re- 
quirement that no individual pion fiber had more than 
3 MeV (89.6%) and the requirement on the target-track 
fitter probability V{xl +xl + X?) (87.4%). 



3. Acceptance factors from 7r SC at events 

Beam pions that scatter in the target had range, en- 
ergy and momentum spectra similar to that of pions from 
K + — > tt + vv and were used to determine the acceptance 
factors associated with the reconstruction and identifi- 
cation of pions in the range stack. Table VI lists the 
acceptance factors measured using 7r scat monitors. 

Candidate events were only accepted if the pion 
stopped in a counter with an operational TD ("Good 
TD" ) . The acceptance factor associated with this require- 
ment was measured on a sample of 7r scat monitor data 
selected by requiring a single pion entering the target, a 
good outgoing track in the target, UTC and range stack, 
a delayed coincidence of less than 5 ns, and range, energy 
and momentum in the signal region. 

In addition to the requirements listed above, 7r scat mon- 
itor data were also required to have a good n — > /i — > e 
signature and the pion was required to stop in a range 
stack counter with an operational TD in order to mea- 
sure the acceptance factor associated with the range stack 
kinematics and tracking. Assignment of target fibers to 
the incoming and outgoing pion in 7r scat was not as robust 
as the assignments made for kaon decays at rest. Mis- 
assignment of target fibers yielded a larger uncertainty 
in the momentum, range and energy calculated for the 
outgoing pion in 7r scat events. The effect of increasing or 
decreasing the signal region in momentum, range and en- 
ergy (Section IIIB4) by ±1 standard deviation was used 
to estimate the systematic uncertainty in this acceptance 
factor. 

The requirements used to assess the acceptance factor 
associated with a non-operational TD in the stopping 
counter were supplemented by requiring a good track in 
the UTC and application of the RNGMOM cut in order 
to measure the acceptance factor associated with the tt — > 
yu — > e signature. The cut on the measured dE/dx in 
range stack counters and the cuts on the consistency of 
the range stack and UTC track were slightly correlated 
with the suite of cuts used to define the 7r — > /i — > e 
signature [55]. The acceptance factor of the tt — > fj, — > e 
signature was assessed both with and without these cuts 
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applied to estimate the systematic uncertainty due to 
these correlations. In addition, a relative correction of 
+1.4% for pion decay-in-flight and pion absorption in the 
stopping counter, estimated from simulation, was applied 
to the acceptance factor associated with the tt — > [i — > e 
signature. 

TABLE VI: Acceptance factors measured with 7r scat monitor 
triggers. The first and second uncertainties are statistical 
and systematic, respectively. The assessment of systematic 
uncertainties is described in the text. A Tscat is the product of 
the three acceptance factors. 

Cut Acceptance factor 

Good TD in stopping counter 0.9984 ± 0.0001 

Range stack kinematics 0.8259 ± 0.0013 ± 0.0120 

tt -> fi -> e signature 0.4805 ± 0.0015 ± 0.0160 

A^ acat 0.3980 ± 0.0014 ± 0.0140 



the same sector in the two innermost layers and in the 
IC. The simulation did not include the acceptance loss 
due to gaps between the neighboring T counters or due 
to insufficient scintillation light in the thin T counters. 
These acceptance losses were measured by using and 
K v 2 decays in KB monitor events. The energy loss in the 
T counter by the charged track differs for and K % 2 
events and simulated events were used to obtain the aver- 
age energy loss for each decay. The measured acceptance 
factors for and were then extrapolated to esti- 
mate 

At-2-ic = 0.9505 ± 0.0012 ± 0.0143 (14) 

where a ±1.5% systematic uncertainty was assigned to 
account for the extrapolation of the UTC track to the T 
counter. 



4- Acceptance factors from simulated events 

Simulated K + — ► tt + vv events were used to evaluate 
the acceptance factors associated with the trigger, phase 
space and 7r + -nuclear interactions in the detector. The 
acceptance of the Ll.l and LI. 2 (DC) components of 
the trigger as described in Section II B were evaluated 
with K V 2 {K^) monitors as described previously in this 
Section. The acceptance of the remaining trigger compo- 
nents is given in the second row of Table VII. The phase 
space acceptance includes the loss due to tt + absorption 
in the stopping counter and decay- in- flight. Aside from 
7T+ absorption in the stopping counter, neither the trig- 
ger nor phase space acceptance include the effect of nu- 
clear interactions. As indicated in Section III A, the com- 
bined trigger and phase space acceptance factor of 11.8% 
(= 0.3225 x 0.3650) was larger than the corresponding 
factor of 6.5% for the -kvv{\) region [2]. The acceptance 
factor associated with 7r + -nuclear interactions was eval- 
uated separately and is given in the fourth row of Ta- 
ble VII. For the -kvv(1) region, the acceptance factor 
due to nuclear interactions was 49.5% [2]. 



6. Normalization to the K^a branching ratio 

We assessed the fraction (f s ) of K + that stopped in the 
target by normalization to the branching ratio [24] 
as described in [2] 

f a = 0.7740 ±0.0011 . (15) 

7. Confirmation of the branching ratio 

The K^2 branching ratio was measured using the K^i 
monitor trigger data in order to confirm the validity of 
the majority of acceptance factors and corrections cal- 
culated with data and simulation. The acceptance fac- 
tors associated with the photon veto (Table IV) were not 
checked by this procedure because the standard photon 
veto cuts were not applied for this measurement. Our 
measurement followed the same analysis procedure as de- 
scribed in [2] but utilized the selection criteria developed 
for the m>P(2) analysis. From this analysis we obtained 



TABLE VII: Acceptance factors determined from simulated 
K+ — > -k + vv decays. Amc is the product of the three accep- 
tance factors. The uncertainties are statistical. 



Component 


Acceptance 


factor 


Trigger 


0.3225 ± 


0.0015 


Phase space 


0.3650 ± 


0.0027 


7r + -nuclear interactions 


0.8284 ± 


0.0104 


Amc 


0.0975 ± 


0.0009 



5. Correction to the T ■ 2 • IC efficiency 

The T ■ 2 • IC component of the trigger (Section II B) 
required a coincidence between range stack counters in 



B(K + -> tt+tt ) = 0.221 ± 0.002 (16) 

where the uncertainty is statistical. The 6% difference 
with the world average value [24] of 0.209 ± 0.001 for the 
branching ratio was taken into account in the assigned 
systematic uncertainty discussed in the next section. 



8. Overall acceptance and sensitivity 

The total acceptance was evaluated as the product 
of A Kjr2 , A Kfi2 , A nscat , A MC , f s and A T .%ic giving 
A tot = (1.37 ± 0.14) x 10~ 3 where we assigned a 10% 
uncertainty on the total acceptance to accommodate the 
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discrepancy in B(K + — > 7r + 7r°) and the additional sys- 
tematic and statistical uncertainties in the acceptance 
evaluated in this Section. Based on the total exposure 
of Nk = 1.71 x 10 12 stopped kaons for this analysis, 
the single event sensitivity (SES = l/(Nj(A tot )) of the 
nvi>(2) analysis was SES = (4.28 ± 0.43) x 1(T 10 which 
can be compared with the SES of the E949 -kvv{1) analy- 
sis of (2.55±0.20)xlCT 10 [2] and the combined SES of the 
previous nvi>(2) analyses of (6.87±0.04) x KT 10 [35, 45]. 



E. Likelihood method 

We determined B{K + — > tt + i>D) using a likelihood 
method that took into account the distributions of the 
predicted background and acceptance within the signal 
region. The signal region was divided into nine cells with 
differing acceptance-to-background ratios as described 
below. The likelihood ratio X was defined as 



XsU e i^L/i-jL (17) 

i=l 

where Sj and bi were the estimated signal and background 
in the z th cell, respectively, di was the observed number 
of signal candidates in the i th cell and n was the total 
number of cells [56]. The estimated signal in each cell 
was given by Sj = B(K + — > 7r + w)/SESi where SESj was 
the single event sensitivity of the i th cell. 

The division of signal region into nine cells was per- 
formed using combinations of the decay pion kinematics 
(KIN), photon veto (PV), delayed coincidence (DC) or 
7T — > \i — > e (TD) cuts. We defined a standard (e.g. 
KINs) and a more restrictive or "tight" version of each 
cut as described previously in Sections IIIB4 (KINt), 
IIIB7 (PV T ), IHB6 (DCt) and HIB5 (TD T ), and 
KINr was defined as KIN S = KIN T +KIN R . The signal 
region was defined by the application of the standard ver- 
sion of all cuts. The signal region was then subdivided 
into cells by the additional selective application of the 
tight version of each cut or the inverted cut as shown in 
Table VIII. 

The additional signal acceptance factor and rejection 
(Table IX) for each of the four tight cuts was deter- 
mined using analogous techniques and samples as de- 
scribed in Sections HID and IIIC, respectively. Based 
on studies of data and simulated events, the background 
components for which no rejection is given in the Table 
were reduced by the acceptance factor of the particular 
cut. For example, the acceptance of the cell defined by 
KINt • TDr • DCr • PVt relative to the acceptance for 
the entire signal region was A KINt x (1 — A TDt ) x (1 — 
4 DCt ) x Ap VT with obvious notation, and the -fT^-target- 
scatter background component relative to the contribu- 
tion to the entire signal region was reduced by the factor 

l/i?KiN T X (1 - A TDt ) X (1 - A DCt ) X l/i?pv T . 



TABLE VIII: The estimated signal-to-background (s/b) and 
background (b) for the nine signal region cells. The cuts 
defining each cell are also given. The s/b is calculated as- 
suming B(K + —> -k + vv) — 1.73 x 1CF 10 . The uncertainties 
on b and s/b are omitted from the table. In the definition 
column, KINs indicates the region defined by the standard 
decay pion kinematics cut and KINt indicates the region de- 
fined by the tight version of the kinematics cut. KINr is 
defined as KINs = KINt + KINr. Analogous designations 
are used for the photon veto (PV), delayed coincidence (DC) 
and -K — > n — > e (TD) cuts. 



Cell 


Definition 


b 


s/b 


1 


KINt 


TD t 


DC T 


PV T 


0.152 


0.84 


2 


KINt 


TDr 


DCt 


PVt 


0.038 


0.78 


3 


KINt 


TDt 


DCr 


PVt 


0.019 


0.66 


4 


KINt 


TDr 


DCr 


PVt 


0.005 


0.57 


5 


KINt 


TD t 


DCt 


PVr 


0.243 


0.47 


6 


KINt 


TDr 


DCt 


PVr 


0.059 


0.45 


7 


KINt 


TD t 


DCr 


PVr 


0.027 


0.42 


8 


KINt 


TDr 


DCr 


PVr 


0.007 


0.35 


9 


KINr 


•TDs 


•DCs 


PVs 


0.379 


0.20 



TABLE IX: Additional acceptance factor (A) or rejection (7?) 
for the tight version of the KINt, PVt, DCt and TDt cuts 
for specific backgrounds. 





KINt 


TD T 


DC T PVt 


A 


0.812 


0.812 


0.911 0.522 


R 


1.63(7^2) 


3.08 (Muon) 


6.3 (CEX) 2.75 (K^) 




2.70(K e4 ) 




1.0 (Beam) 2.75 {K„ 2 y) 




1.20(^2 7 ) 







IV. RESULTS 

In this Section, we describe the results of exam- 
ining the signal region and the evaluation of the 
K + — ► tt + vv branching ratio. We also describe the eval- 
uation of our observations within alternative models of 
K + — ► 7T++ nothing. 



A. Examination of the signal region 

After completion of the background and acceptance 
analyses, all selection criteria were applied to the 
•kvv{1+2) trigger data. Three signal candidate events 
were observed. Some measured properties of the three 
events are listed in Table X. The plot of the kinetic en- 
ergy vs. range of the three events along with the events 
found in the previous 7iw(l) [2] and ttvD(2) [35, 45] anal- 
yses are shown in Figure 16. The candidates' measured 
properties used in the selection criteria were consistent 
with the expected distributions for K + — > tt + vv decay 
signal. There was no observed activity in the kaon fibers 
at the time of the 7r + for any of the three candidates ac- 
cording to the CCDPUL analysis. All three events failed 
one or more of the tight cuts described in Section HIE. 
The 7r + momentum of event A failed the tight momen- 
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TABLE X: The cell (Table VIII) and the momentum, range 
and energy of the 7r + in the three signal candidate events A, 
B and C. The measured decay times are in units of K + , tt + 
and lifetimes [24], appropriately. 



Event 


A 


B 


C 


Cell 


9 


5 


7 


Momentum (MeV/c) 


161.5 


188.4 


191.3 


Range (cm) 


17.3 


24.2 


26.1 


Kinetic energy (McV) 


76.1 


95.6 


97.9 


K + decay time 


0.30 


1.27 


0.42 


7r + decay time 


0.86 


0.64 


0.39 


H + decay time 


2.71 


1.03 


4.33 



turn cut of 165 MeV/c. Events B and C failed the tight 
photon veto cut due to energy deposits of 2.4 and 2.1 
MeV in the EC above the 1.7 MeV threshold, respec- 
tively, and events A and C failed the tight DCt cut of 6 
ns on the kaon decay time. 



The K+ -> 7T+ vv branching ratio 



+ vv branching ratio 
-> it + vv) that max- 



The central value of the K + — » tt 
was taken to be the value of B{K + 
imized the likelihood ratio X (Equation (17)). For the 
three events observed by this analysis, we determined 
B(K+ -> tt+vv) = (7.89±§;^) x lfT 10 where the quoted 
68% confidence level interval was determined from the 
behavior of X as described in [56] and took into account 
both the statistical and systematic uncertainties. The 
systematic uncertainties included the 10% uncertainty in 
the acceptance as well as the uncertainties in the esti- 
mation of the background components. The inclusion 
of systematic uncertainties had a negligible effect on the 
confidence level interval due to the poor statistical pre- 
cision inherent in a three event sample. The probabil- 
ity that these three events were due to background only, 
given the estimated background in each cell (Table VIII) , 
was 0.037. 

When the results of the previous -nvviV) and 7rz/p(2) 
analyses [2, 35, 45] were combined with the results 
of this analysis, we determined B(K + — > tt + vv) — 
(1. 73^05) x 10~ 10 [57]. Systematic uncertainties were 
treated as described above when performing the combi- 
nation, except that we assumed a correlated 10% uncer- 
tainty for the acceptance assessed by each analysis. The 
probability that all seven events were due to background 
only (background and SM signal) was estimated to be 
0.001 (0.073). 

The K + — > -k + vv branching ratio for each of two pion 
momentum regions are given in Table XI along with the 
SM predictions. The combined E787 and E949 wvv(2) 
data were used for the branching ratio in the 
[130,205] ([205,227]) MeV/c region. The limits of the 
two momentum regions were determined by the require- 
ments on the reconstructed pion region, the detector reso- 
lution and the desire to have contiguous, non-overlapping 
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FIG. 16: Kinetic energy vs. range of all events passing all 
other cuts. The squares represent the events observed by this 
analysis. The circles and upward-pointing triangles repre- 
sent the events observed by the E787 and E949 nvv{l) anal- 
yses, respectively. The downward-pointing triangle represent 
the events observed by the E787 nvv(2) analyses. The solid 
(dashed) lines represent the limits of the nvv(l) and -kuv{2) 
signal regions for the E949 (E787) analyses. Despite the 
smaller signal region in R w vs. E^, the ixvviV) analyses were 
4.2 times more sensitive than the nvi>(2) analyses. The points 
near E w = 108 MeV were decays that survived the pho- 
ton veto cuts and were predominantly from the m/i>(l) analy- 
ses due to the higher sensitivity and the less stringent photon 
veto cuts. The light gray points are simulated K + — * n + ui? 
events that were accepted by the nvv(l+2) trigger. 



regions. The boundary between the two regions of 205 
MeV/c was determined by the lower and upper limits 
on the reconstructed ir + momentum that were set to be 
approximately 2.5 standard deviations from the nominal 
K„2 momentum for the Truu(l) and 7wP(2) analyses, re- 
spectively. The lower limit of 130 MeV/c was determined 
by the irvD(2) requirement on the reconstructed pion mo- 
mentum. The upper limit of 227 MeV/c is the kinematic 
limit for the K + — > it^vv decay. 

The 90% CL upper limit B(K+ -> ir+vv) < 
3.35 x 10~ 10 was also determined and can be used to cal- 
culate a model-independent upper limit of 14.6 x 10~ 10 
on the CP- violating process K® — > 7r°i/p[58]. This limit 
is substantially smaller than the current experimental 



limit of B(K9 



< 670 x 10 
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TABLE XI: The measured branching ratios in units of 1CP 10 
assuming a standard model, scalar or tensor interaction [57] 
in two momentum regions as described in the text. The SM 
prediction was taken from [5] and scaled to the two pion mo- 
mentum regions using the SM spectral shape for K + — > -k + vv. 
The 68% and 90% CL intervals for the double ratio, described 
in the text, are also given. 



Momentum 
range (MeV/c) 


SM 
prediction 


Interaction 
SM Scalar Tensor 


[130, 205] 
[205, 227] 


0.49 ± 0.04 
0.28 ± 0.02 


2.9ltt79 
0.49t°;« 


0.591^ 
3.50t!:» 


0.49+— 
3.103$ 


Double ratio 68% CL 
90% CL 


(0.56,6.15) 
(0.24,17.6) 


(0.29,3.32) 
(0.13,9.55) 


(0.08,0.88) 
(0.04,2.66) 



C. Alternative model interpretations 

The combined results of the E787 and E949 experi- 
ments can also be interpreted assuming a scalar or ten- 
sor interaction. The spectra are compared with the 
SM spectrum in Figure 17. Using the same treatment 
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FIG. 17: (Top) The generated 7r + kinetic energy spectrum for 
SM (solid), scalar (dashed) and tensor (dotted) interactions. 
(Bottom) The distributions for events passing the -kvv(1+2) 
trigger. 

of the data as for the standard model interpretation, 
the branching ratio for the scalar and tensor spectra 
alone were B SC (K+ -> tt+vv) = (9.9±|;|) x 1(T 10 and 
B tc (K+ -> n+vv) = (4.9±|;|) x lO" 10 , respectively, or 
B SC (K+ -> tt+vv) < 21 x 1CT 10 and B tc (K+ -» n+vv) < 
10 x 1CT 10 at 90% CL. 

The consistency of the distribution of the observed 
events with the shape of the SM, scalar or tensor spec- 



trum was evaluated using the double ratio 

£(130,205) / £(205,227) 
~ JT(130, 205) / ^(205, 227) 1 ' 

where B(P\,P2) is the branching ratio and J r (Pi,P 2 ) 
is the fraction of phase space in the momentum range 
(Pi,P 2 ) MeV/c. The ratio .F(205, 227)/J"(130, 250) is 
0.084, 0.014 and 0.59 for the scalar, tensor and stan- 
dard model spectra, respectively. If the distribution of 
observed events were consistent with the shape of the as- 
sumed spectrum, then the double ratio would be unity. 
The 68% and 90% CL intervals of the DR and the branch- 
ing ratios are given in Table XL 

The data have also been interpreted in the two-body 
decay model, K + — ► 7T+AT, where A" is a massive non- 
interacting particle, either stable or unstable. The 90% 
CL upper limit on the branching ratio as a function of 
the mass of X is shown in Figure 18. For the case of an 
unstable AT, the decay of X was assumed to be detected 
and vetoed with 100% efficiency if X decayed within the 
outer radius of the BV. The E949 limit of B{tt° -> vv) < 
2.7 x 10~ 7 at 90% CL [60] can be combined with the 
world average value of B(K + — > 7r+7r°) [24] to set a 90% 
CL limit of B(K+ -> tt+AT) < 5.6 x 10~ 8 for M x = M n o 
with X stable. 




Mass (MeV/c ) 

FIG. 18: The solid line represents the 90% CL upper limit on 
B{K + — > ir + X) as a function of the mass of X assuming X 
is stable. The dashed lines represent the 90% CL for X with 
the indicated lifetimes as described in the text. 

The limit in the Figure 18 can also be interpreted as 
a limit on the product of branching ratios B{K + — > 
7T + P°) x B(P° — > vv) for a hypothetical, short-lived 
particle P°. The HyperCP collaboration observed three 
events consistent with E + — > pP° with P° — > hav- 
ing a mass M(P°) = 214.3 ± 0.5 MeV/c 2 [61]. A mass 
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of 214.3 MeV/c 2 would correspond to a recoiling n + mo- 
mentum, range and energy of 170.1 MeV/c, 19.5 cm and 
80.5 MeV, respectively, in a two-body K + decay. Of the 
four events in the -kvv{2) region observed in E787 [45] and 
E949, the closest was the present candidate A that dif- 
fered from the expected 7r + momentum, range and energy 
by 3.7, 2.4 and 1.5 standard deviations, respectively. 



V. CONCLUSION 
A. Summary 

The branching ratio of the decay K + — > ir + vv is 
precisely predicted in the standard model to be 
B(K+ -» tt+vu) = (0.85 ± 0.07) x lO" 10 [5]. As the de- 
cay is a flavor-changing neutral current process and thus 
sensitive to new physics effects [6], it represents an op- 
portunity to unambiguously probe for physics beyond the 
standard model. 

The E787 and E949 experimental program at 
Brookhaven National Laboratory has demonstrated the 
feasibility of observation of the rare decay K + — > tt + vv 
using stopped kaons despite the challenging experimen- 
tal signature. The analysis presented in this article has 
established that backgrounds can be reduced to a rea- 
sonable level while maintaining signal acceptance to en- 
able the irvv(2) region to be a viable supplement to the 
ttvv(\) region [2] in the measurement of the K + — ► ir + vv 
branching ratio. The branching ratio B(K + — ► ir + vv) — 
(1.73±i jf) x lO" 10 for the seven events observed by E787 
and E949 is consistent with the SM expectation. 



B. Final comments 

We briefly summarize some important issues that 
emerged from the E787 and E949 searches for 
K + — » n + vD and the measurement of B(K + — > 7r + ^P) in 
a stopped K + experiment: 

1. Importance of blind analysis. The analysis proce- 
dure of concealing or obscuring the contents of the 
signal region [33, 62] are now well-established and 
widespread in particle physics. The main benefit is 
to avoid or minimize bias in the selection criteria. 
This is particularly important for a rare process 
such as K + — > tt + vv with a poor experimental sig- 
nature that requires many individual cuts. 

2. Use of data to estimate background and accep- 
tance. In conjunction with the blinding of the sig- 
nal region, the division of the data into 1/3 and 2/3 
samples and the use of two powerful independent 
cuts provides background estimates that take into 
account instrumental effects not present in simula- 
tion and provides sensitivity to background beyond 
that of the K + — ► n + vv sample as first noted in 



[33] . The bifurcation technique also permits exam- 
ination of events that occur near the signal region 
to validate the background estimates and investi- 
gate unforeseen sources of background. 

3. Unforeseen acceptance losses. In the E787 pro- 
posal [63], only the ixvv(X) analysis region was con- 
sidered to be viable and the estimated acceptance 
was 1.5%. This can be compared to the E949 
-nvv(V) acceptance of 0.22% [2] and 0.14% for the 
present analysis (Section HID). Accidental activ- 
ity reduced the acceptance factors associated with 
vetoing. Acceptance was also reduced due to the 
cuts to suppress background, particularly related 
to muons in the irvv(l) region and pion scatter- 
ing in the irvv(2) region, that needed to be more 
stringent than anticipated as well as due to cuts to 
suppress beam-related background that were not 
considered in the proposal. 

4. Importance of redundancy in detector systems. In 
E949 and its predecessor E787, since nearly every 
detector element participated in the veto function 
for dealing with additional particles, the loss of any 
element represented a reduction in background re- 
jection and motivated the redundant use of ADCs, 
TDCs and high-speed waveform digitizers (CCDs 
or TDs). 

5. Need for 4-7T sr photon veto coverage. The early 
K + — > ir + v9 counter experiments showed the ben- 
efits of photon veto capability over the full 4-7r sr 
solid angle [27, 28]. Each modification or upgrade 
of the original E787 experiment included photon 
veto enhancements. For example, E949 contained 
upgrades with respect to E787 that sought to fur- 
ther suppress the contribution of the K^i back- 
ground to both the irvv(l) and ttvv(2) analyses. 
The barrel veto liner (BVL) improved the photon 
veto rejection in the barrel region by a factor of 2 
and the active degrader (AD) increased the photon 
veto rejection for the ttvv(2) analysis by 2 (Sec- 
tion IIIB7). 

Extrapolating from the E787 and E949 experience, it 
would be possible to extend the stopped K + experi- 
ment approach to yield orders of magnitude more ob- 
served K + — ► ir + vv decays [64]. As indicated above, the 
high rate of kaon interactions at 710 MeV/c resulted in 
reduced acceptance due to accidental spoiling of good 
events and increased backgrounds requiring highly re- 
strictive cuts. If a higher primary proton beam inten- 
sity were available (with a high duty factor), lower mo- 
mentum kaons, e.g. 450 MeV/c, would result in a much 
higher stopping fraction and reduced interactions which 
cause accidentals. To take advantage of this effect, a 
shorter particle-separated beam line (e.g. 13 m compared 
to the present 19 m), would be necessary to reduce atten- 
uation of the kaon beam due to decay. Other improve- 
ments which would result in increased acceptance might 
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involve use of a higher magnetic field (IT — ► 3T) allow- 
ing a more compact detector with improved momentum 
resolution, finer segmentation of the range stack, a more 
hermetic, thicker fully active veto detector (e.g. crystals 
or liquid xenon), and an improved target. 
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